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PREFACE 


This  technical  report  consists  of  publications  and  preprints  which  are 
representative  of  research  supported  by  this  contract  from  September  1982  thru 
August  1983.  The  purpose  of  this  preface  is  to  comment  on  the  motivation  of, 
and  background  for,  this  research  and  to  interrelate  the  topics  covered. 

Paper  No.  1  describes  the  first  detailed  observations  of,  and  theory 
for,  focused  acoustical  backscattering  from  an  elastic  sphere.  This  focusing 
is  referred  to  as  "axial  focusing"  because  the  scattered  pressure  is  locally 
enhanced  near  the  backward  axis  as  evidenced  by  Fig.  9-12.  This  phenomenon  is 
called  the  "glory"  since  the  enhanced  optical  backscattering  from  cloud 
droplets  (also  known  as  the  glory)  is  caused  in  part  by  axial  focusing.  The 
first  detailed  model  of  transmitted  wave  glory  amplitudes  from  spheres 
(acoustical  or  optical)  was  developed  under  this  contract.1  Paper  No.  1 
contains  several  theoretical  extensions  of  our  previous  (fluid  sphere)  model1 
and  it  gives  data  which  confirm  those  extensions.  Axially-focused 
backscattering  is  potentially  significant  to  sonar  technology.  Due  to  the 
Intrinsic  isotropy  of  spheres,  there  has  long  been  an  interest  in  the  use  of 
fluid-filled  spherical  shells  as  sonar  calibration  targets.  These  targets  may 
also  be  useful  as  passive  navigational  beacons  for  which  backscattering  is 
large  for  a  broad  range  of  frequencies.  The  theoretical  considerations 
described  should  facilitate,  for  the  first  time,  a  detailed  understanding  of 
focused  echoes  from  large  spheres  having  ka  »  1. 

Paper  No.  2  describes  a  detailed  calculation  of  the  temporal 
signatures  of  axially-focused  backscattering  from  spheres.  It  is  shown  that 
the  time-domain  response  due  to  axial  focusing  is  descrlbable  (under 
conditions  noted)  by  the  "half-order  derivative"  operator.  The  effect  of  this 
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operator  on  a  sine-wave  tone  burst  is  calculated  and  discussed.  The  quasi- 
steady-state  approximation  used  in  Paper  No.  1  is  shown  to  be  applicable  for 
the  conditions  of  that  experiment.  In  addition,  axially- focused 
backscattering  from  a  large  torus  is  modeled.  The  frequency  domain  response 
associated  with  axial  focusing  appears  to  be  one  of  the  canonical  response 
functions  of  acoustics  so  the  temporal  response  is  applicable  to  several  other 
problems  in  acoustics,  some  of  which  are  noted  below. 

Paper  No.  2  also  reviews  Lighthill's  formulation  of  the  pressure 
radiated  from  a  line  source  in  terms  of  the  half-order  derivative  operator.  A 
specific  example,  the  transient  radiation  from  certain  optoacoustic  sources, 
is  considered.  When  the  optical  intensity  modulation  i3  that  of  a  half-cycle 
of  a  sine  wave,  the  calculation  displayed  in  Fig.  5  should  approximate  the 
radiated  pressure.  This  calculation  (and  similar  ones  for  other  inputs) 
should  be  useful  for  calibrating  broad-band  hydrophone  systems  needed  in  a 
study  (supported  by  this  contract)  of  other  optoacoustics  sources. 

Paper  No.  3  describes  the  first  uniform  analysis  (i.e.  an  analysis 
which  is  equally  valid  for  all  scattering  angles)  of  a  novel  measurement 
configuration  for  light  scattering.  The  measurement  configuration  considered 
is  potentially  useful  for  the  discrimination  of  light  scattered  from  spheres 
(small  air  bubbles  in  water  or  fog  droplets)  from  that  of  non-spherical 
particulates.  It  was  previously  introduced  by  us  to  facilitate  the  detection 

of  optical  glories  (i.e.  axialy-focused  backscattering)  from  air  bubbles  in 

2  3  '  n  s 

liquids.  ’  Our  previous  computations  *  and  unpublished  visual  observations 
indicate  that  optical  glories  enhance  the  near-backward  scattering  from  small 
air  bubbles  freely  rising  in  water.  (Small  bubbles  significantly  affect  the 
acoustical  properties  of  water  and  also  the  air-sea  interaction  at  the  ocean's 
surface.)  The  calculation  described  is  applicable  to  the  production  of  cross- 
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polarized  light  due  to  scattering  from  small  bubbles  (and  from  fog  or  cloud 

droplets)  provided  the  number  density  of  scatterers  is  sufficiently  small  that 

multiple  scattering  may  be  neglected.  The  calculation  is  also  applicable  to 

non-cross-polarized  (sometimes  known  as  copolarized)  scattering  and  to  forward 

as  well  as  backward  optical  glories. 

In  other  work  supported  by  this  contract  we  have  observed  the 

6  7 

scattering  of  both  white  light0*  and  laser  light  from  bubbles  in  the  vicinity 

of  the  critical  scattering  angle.  Color  photographs  of  the  white-light 
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scattering  are  to  be  published. 

In  preparing  the  papers  enclosed  in  this  report  the  principal 
investigator  (P.  L.  Marston)  has  benefited  from  the  collaboration  with,  and/or 
critical  comments  by,  several  students.  These  include:  T.  J.  B.  Hanson,  D. 

S.  Langley,  B.  T.  Unger,  and  K.  L.  Williams. 


P.  L.  Marston 
September  1983 
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Paper  Number  1: 

Observation  of  the  acoustic  glory: 

High-frequency  backs cat te ring  from  an  elastic  sphere 

a) 

Philip  L.  Mar s ton,  Kevin  L.  Williams,  and  Timothy  J.  B.  Hanson 

Department  of  Physics 
Washington  State  University 
Pullman ,  Washington  99164 

Reprinted  from:  J.  Acoust.  Soc.  Am.  74.  605-618  (1983). 

Measurements  of  the  scattering  of  2.10  MHz  ultrasonic  bursts  gave  direct 

evidence  of  the  acoustic  glory.  Echo  amplitudes  were  measured  as  a  function 

of  the  transverse  displacement  x  from  the  backward  axis;  they  were  maximized 

on  the  axis  and  they  manifested  side  lobes.  Except  near  minima,  the  amplitudes 

were  proportional  to  J_(kbx/z  )  where  b  and  z  are  the  radius  of,  and  axial 

u  n  n 

distance  to,  the  virtual  focal  circles  from  which  the  echoes  appeared  to  emanate. 
This  dependence  on  x  is  direct  evidence  of  axial  focusing  (The  principal 
feature  of  optical  glories  of  drops  and  bubbles  is  also  axial  focusing.).  The 
experiment  was  performed  in  water  with  a  fused  silica  sphere  having  a  radius 
a=51.6  mm  for  which  2Tra/X  ■  ka=457.  A  theory  is  developed  for  certain  of  the 
echo  amplitudes  and  focal  parameters.  The  theory  typically  overestimates  the 
amplitudes  of  the  stronger  of  the  echoes  by  5Z.  Appendices  describe  aspects 
of  the  theory  not  amenable  to  the  experiment,  including:  (i)  a  shift  of 
focal  parameters  which  should  be  significant  when  ka£l00,  and  (ii)  a  distortion 
of  transients  described  in  part  by  the  half-order  derivative  operator. 

a)  Present  address:  Defense  Systems  Division, 

Honeywell,  Inc.,  Hopkins,  MN  55343 

PACS  Nos:  43.20. Fn,  43.20.Px,  43.20.Dk,  42. 10. He 
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Observation  of  the  acoustic  glory:  High-frequency  backscattering 
from  an  elastic  sphere 

Philip  L.  Marston,  Kevin  L.  Williams,  and  Timothy  J.  B.  Hanson** 

Department  of  Physics,  Washington  State  University,  Pullman,  Washington  99164 

(Received  28  December  1982;  accepted  for  publication  9  May  1983) 

Measurements  of  the  scattering  of  2. 10-MHz  ultrasonic  bursts  gave  direct  evidence  of  the  acoustic 
glory.  Echo  amplitudes  were  measured  as  a  function  of  the  transverse  displacement  x  from  the 
backward  axis;  they  were  maximized  on  the  axis  and  they  manifested  side  lobes.  Except  near 
minima,  the  amplitudes  were  proportional  to  J0(kbx/z„ ),  where  b  and  z„  are  the  radius  of,  and 
axial  distance  to,  the  virtual  focal  circles  from  which  the  echoes  appeared  to  emanate.  This 
dependence  on  x  is  direct  evidence  of  axial  focusing.  (The  principal  feature  of  optical  glories  of 
drops  and  bubbles  is  also  axial  focusing.)  The  experiment  was  performed  in  water  with  a  fused 
silica  sphere  having  a  radius  a~5 1.6  mm  for  which  2 ira/A  —  £a~457.  A  theory  is  developed  for 
certain  of  the  echo  amplitudes  and  focal  parameters.  The  theory  typically  overestimates  the 
amplitudes  of  the  stronger  of  the  echoes  by  5%.  Appendices  describe  aspects  of  the  theory  not 
amenable  to  the  experiment,  including:  (i)  a  shift  of  focal  parameters  which  should  be  significant 
when  ka  $  100,  and  (ii)  a  distortion  of  transients  described  in  part  by  the  half-order  derivative 
operator. 

PACS  numbers:  43. 20.Fn,  43.40.Ey,  43.20.Px,  42.10.Hc 


INTRODUCTION 

There  have  been  numerous  investigations  of  acoustic 
backscattering  from  solid  elastic  spheres  in  water. 1-10  (Much 
of  this  work  is  reviewed  in  Ref.  8.)  When  the  sphere  is  very 
large  in  comparison  to  an  acoustic  wavelength,  a  simple  de¬ 
scription  of  the  backscattering  should  be  possible  by  ac¬ 
counting  for  reflection  and  the  various  longitudinal  and 
shear  rays  internal  to  the  sphere.*"'2  It  is  insufficient  to  use 
elementary  ray-optical  methods10  because  of  a  novel  focus¬ 
ing  effect. 13  Diffraction  must  be  included913  in  models  of  this 
“axial  focusing.”  It  is  the  purpose  of  this  paper  to  describe 
observations  of  axial  focusing  of  backscattering  and  to  ex¬ 
tend  a  previous  model,13  developed  for  fluid  spheres,  to  the 
case  of  the  backscattering  of  transient  signals  from  elastic 
spheres.  It  is  demonstrated  that  the  scattering  to  angles 
slightly  displaced  from  the  exact  backward  direction  can  be 
greatly  influenced  by  the  size  of  the  scatterer. 

It  is  appropriate  to  refer  to  the  scattering  described'here 
as  the  acoustic  “glory”  because  the  term  “glory”  is  generally 
applicable  to  backscattering  from  spheres  (or  other  spheri¬ 
cally  symmetric  potentials)  when  enhanced  by  axiai  focus¬ 
ing.  13-13  Axial  focusing  enhances  the  amplitudes  due  to  rays 
backscattered  with  nonzero  impact  parameters  (the  “glory 
rays”)  but  it  does  not  affect  the  amplitudes  of  the  axial  reflec¬ 
tions. 13  There  is  no  axial  focusing  of  the  backscattering  from 
cylinders9-13  and  consequently  the  exact  scattering  (normal¬ 
ized  to  that  of  the  geometric  reflection  firom  a  fixed-rigid 
object  of  the  same  size)  is  somewhat  greater  for  large  elastic 
spheres  than  for  cylinders.6,9 

Figure  1  gives  a  conceptual  description  of  the  experi¬ 
ment.  A  pistonlike  transducer  of  diameter  d  =  12.7  mm  is 
driven  electrically  so  as  to  produce  a  tone  burst  of  short  du- 


*' Proem  address:  Defense  Systems  Division,  Honeywell,  Inc.,  Hopkins, 
MN  33343. 


ration.  The  principal  burst  (neglecting  buildup  and  decay  of 
the  transducer’s  response)  is  four  cycles  of  a  2. 1-MHz  sine 
wave.  This  corresponds  to  a  wavelength  A  of  0.71  mm.  The 
distance  zp  from  the  transducer  to  the  near  point  of  the 
sphere  is  171  cm  which  is  somewhat  greater  than  id2/4A; 
consequently,  the  sphere  is  in  the  farfield  of  the  transducer. 
The  sphere  is  synthetic  fused  silica16  with  a  radius 
a  =  51.363  mm.  It  subtends  a  half-angle  r^a/zp  radian  (rela¬ 
tive  to  the  transducer)  which  is  less  than  the  angle  of  the  first 
minimum  (~  1 .22  A /d  rad)  of  the  farfield  radiation.  The  scat¬ 
tering  in  near  backward  directions  was  monitored  with  a 
hydrophone  (at  the  end  of  a  bent  hypodermic  needle).  The 
hydrophone  could  be  scanned  along  a  line  transverse  to  the 
symmetry  axis.  The  on-axis  distance  z  to  the  sphere  was  al¬ 
ways  <ka2  where  k  =  2 ir/A.  Consequently,  the  hydrophone 
was  not  in  the  farfield  scattering  region  (in  the  terminology 
of  Ref.  13).  For  this  experiment  fa»~457.  Typically  z~74 
cm  and  displacements  of  the  hydrophone  from  the  axis  by  an 
amount  x  of  a  few  millimeters  significantly  reduced  the  glory 
scattering  amplitudes. 

Principal  modifications  of  the  fluid-sphere  model  (Ref. 
13)  needed  for  the  description  of  the  experimental  results  are 


FIG.  I  Simplified  dUfrem  of  the  acattermf  experiment.  The  hydrophone 
mij  tit  ii  ■mini  »lmi|  i  llnr  mnenrw  in  ll»  qisninr  site  ihflnwl  trr  thr 
source  and  the  qihere.  Tbe  diagram  to  not  drawn  to  scale. 
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description  of  rays  and  focal-circle  parameters  for  an  elastic 
sphere;  use  of  reflection  and  transmission  coefficients  for  wa¬ 
ter-quartz  interface;  modeling  of  nearfield  diffractive  ef¬ 
fects;  and  discussion  of  the  transient  response.  The  modifica¬ 
tions  are  first  developed  for  the  steady-state  response.  Ray 
acoustics  is  used  to  model  amplitudes  at  planes  near  the 
sphere.  These  amplitudes  then  propagate  (with  diffraction) 
to  the  hydrophone.  Equation  numbers  from  Ref.  13  will  be 
cited  with  an  F  prefix.  Our  nomenclature  extends  that  of 
Ref.  13. 

l.  RAYS  AND  GEOMETRIC  FOCAL-CIRCLE 
PARAMETERS 

Backscattered  ray  paths  through  large  elastic  cylin¬ 
ders"12  and  spheres10  depend  on  the  following  integer  pa¬ 
rameters:  n,  the  total  number  of  ray  traverses  of  the  interior; 

m,  the  number  of  traverses  as  a  shear  wave;  n-m,  the  number 
of  traverses  as  a  longitudinal  wave;  and  /,  where  2/— 1  is  the 
number  of  times  internal  chords  cross  the  axis  shown  in  Fig. 
1.  (The  parameter  ri  of  Ref.  13  is  /  —  1.)  To  simplify  the 
description  of  focal-circle  parameters,  the  following  ap¬ 
proximations  will  be  made:  (i)  the  incident  wave  will  be  taken 
to  be  a  plane  wave,  and  (ii)  the  geometric  condition  on  glory 
rays  will  be  that  they  leave  the  sphere  traveling  parallel  to  the 
axis.  Condition  (i)  is  an  approximation  because  zp/a  is  finite 
while  (ii)  is  an  approximation  because  z/a  is  also  finite.  These 
approximations  facilitate  the  use  of  geometric  parameters 
from  Ref.  1 3 .  As  in  Ref.  1 3  the  axial  rays  which  contribute  to 
the  scattering  along  the  backward  axis  have  a  vanishing  im¬ 
pact  parameter  while  the  impact  parameters  are  finite  for 
glory  rays. 

When  the  ray  traverses  are  either  all  longitudinal 
(m  —  0)  or  all  shear  [m  =  /?),  they  will  be  referred  to  as  “un¬ 
mixed”  chords.  Geometric  parameters  are  as  described  in 
Ref.  13  with  the  acoustic  refractive  index  M  replaced  by 
ML  =  cw/cL  and  Ms  =  cw/cs  for  the  respective  cases 
where  cL  and  cs  are  the  bulk  longitudinal  and  shear  velocities 
for  fused  silica  and  cw  is  the  sound  velocity  in  water.  These 
indices  depend  slightly  on  temperature  (see  Sec.  VII);  a  typi¬ 
cal  temperature  was  22  'C  which  had  ML  =0.2511  and 
Ms  =  0.3957.  The  form  of  ray  paths  with  /  =  1  is  illustrated 
in  Fig.  2  which  is  drawn  for  M  =  0.6.  Figure  2  also  serves  to 
define  several  angles  and  geometric  parameters  derived  in 
Ref.  13.  To  an  external  observer  located  near  the  backward 
axis,  the  scattering  appears  to  emanate  from  a  virtual  ring¬ 
like  source,  the  focal  circle. 1 1  The  circle  has  a  radius  b  and  is  a 
distance  a  behind  the  dashed  vertical  plane,  the  exit  plane. 
(The  subscript  n  of  Ref.  13  will  be  omitted;  here  a  and  b 
depend  on  n,  m,  and  /.)  The  angle  of  incidence  6  which  satis¬ 
fies  the  glory  condition  of  exact  backscattering  will  be  de¬ 
noted  by  0„„  i .  For  the  unmixed  cases,  it  is  given  by  Eq. 
,(F17)when/=  1  and  n  =  3  or  4  and  by  the  numerical  proce¬ 
dure  described  in  Ref.  13  when  n  >  4.  The  focal-circle  radius 
is  b  —  ba  where 

b  =  sin  (1) 

a  is  given  by  Eq.  (FlOb).  Several  unmixed  chords,  for  the  case 
of  fused  silica  in  water,  are  shown  in  Fig.  3. 

When  the  chords  are  mixed,  the  glory  condition  on  the 
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FIG.  2.  Backscattered  rays  in  a  sphere  representative  of  unmixed  chord 
glories.  When  the  figure  is  rotated  about  the  CC '  axis,  the  Fn  trace  out  focal 
circles  which  are  the  apparent  sources  of  backward  propagating  toroidal 
wave  fronts.  These  waves  are  axially  focused  along  the  backward  axis.  The 
specular  reflection  appears  to  diverge  from  A„  and  is  unfocused. 


angle  of  incidence  6  is  again  found  by  considering  rays  trav¬ 
eling  parallel  to  the  axis  when  they  are  incident  on  the  sphere 
and  exiting  from  it.  As  in  the  unmixed  case,  b  describes  all  of 
the  following:  the  radius  of  the  virtual  source  (the  focal  cir¬ 
cle)  and  the  ray’s  distance  from  the  axis  at  impact  and  when 
in  the  exit  plane.  The  generalized  Snell’s  law  and  the  spheri¬ 
cal  geometry  give  the  following  condition10'12  on  b 

ir[l  —  j)  =  (n  —  m)arccos(6  /M, ) 

+  m  arccosj/t  /Ms )  -  arccos(6 ).  (2) 

This  has  been  solved  numerically  and  the  glory  ray’s  angle  of 


FIG.  3.  Representative  unmixed -chord  glory  rays  for  a  fuaed-sihca  sphere 
in  water.  The  parameters  (n,  m,  / )  are  indicated.  Dashed  and  solid  rays  rep¬ 
resent  shear  and  longitudinal  waves,  respectively. 
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FIG.  4.  Representative  mixed-chord  glory  rays  for  a  fused  silica  sphere  in 
water.  The  upper  and  lower  groupsshow  all  sequences  forthe(3,  1. 1  (and  (4, 
2,  1  leases,  respectively. 

incidence  found  from  ( 1 ).  Snell’s  law  gives  the  following  con¬ 
dition  on  the  internal  angles  v,  and  vs  (relative  to  the  normal 
at  each  vertex,  see  v  in  Fig.  2)  for  the  longitudinal  and  shear 
chords,  respectively, 

Ml  sin  vL  =  Ms  sin  vs  =  b.  (3) 

There  are  n\/[m\[n  —  m|!]  different  ways  of  ordering  the 
chords. "  u  Figure  4  illustrates  all  choices  for  the  cases  (n,  m, 
/)  =  ( 3,  1,  1 )  and  (4,  2,  1). 


Consider  the  propagation  phase  shift  between  the  cross¬ 
ing  of  the  dashed  vertical  plane  (Fig.  2)  of  an  incoming  ray 
and  the  second  crossing  of  that  plane  as  the  ray  exits.  For  the 
exactly  backscattered  ray,  this  phase  shift  is 

=2ka[\  —  cos  6„m, 

+  mMs  cos  vs  +  (n  -  m]MLco&  vL  ]  (4) 

in  both  the  mixed  and  unmixed  cases.  For  the  mixed  cases, 
Vn.mj  (and  the  associated  pulse  delay)  are  independent  of  the 
ordering  of  the  shear  and  longitudinal  chords.1 112  Section  II 
demonstrates,  however,  the  combined  reflection-transmis¬ 
sion  coefficient,  needed  for  the  description  of  scattering  am¬ 
plitudes,  depends  on  the  ordering. 

A  dimensionless  parameter  q,  describing  the  spreading 
of  the  wave  front  (9  is  the  limit  of  the  ratio  of  the  arc  lengths 
of  d  'e'  and  de  in  Fig.  2  as  the  latter  vanishes),  is  needed  to 
model  the  scattering.13  For  unmixed  chords,  q  is  given  by 
Eq.  (F12).  For  the  mixed  chord  cases,  neither  a  y:  q  have 
been  computed  due  to  the  complexities  required  to  extend 
the  derivation  given  in  Appendix  A  of  Ref.  13.  Consequent¬ 
ly,  the  description  of  mixed-chord  scattering  will  be  limited 
to  modeling  of  the  dependences  on  scattering  angle  and  re¬ 
flection-transmission  factors. 

II.  REFLECTION  AND  TRANSMISSION  FACTORS 

To  model  the  strength  of  the  virtual  ringlike  source  at 
each  focal  circle,  the  procedure  from  Ref.  13  will  be  fol¬ 
lowed.  Reflection  and  transmission  coefficients  appropriate 
to  an  unbounded  plane  surface  are  used  in  conjunction  with 
geometric  factors  and  energy  arguments.  Consider  an  inci¬ 
dent  traveling  plane  wave  carrying  an  energy  Ein  in  a  speci¬ 
fied  time  interval.  After  a  reflection  from,  or  transmission 
through,  the  interface  let  the  corresponding  temporal  por- 


FIG.  5.  Square  root  of  the  energy  ratios 
for  reflection  and  transmission-fol- 
lowed-by-reflection  These  were  com¬ 
puted  for  a  plane  water-fused  silica  in¬ 
terface.  The  ray  in  the  water  (which  gives 
rise  to  L  and  S internal  rays)  has  an  angle 
of  incidence  8. 


W7 


8  (degrees) 
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TABLE  I.  Reflection  and  focal  parameters  and  form  functions  for  backscattering  for  selected  unmixed-chord  glories.  The  sphere  is  filled  silica  in  fresh  water 
at  22*C. 


(a,  m,  /) 

^«./(4eg) 

IBi 

Tw,R, 

B 

a/a 

9 

b /a 

|C,  (457, 0|| 

3,0,1 

7.77 

0.27 

0.14 

2.1£— 2 

1.038 

27.3 

0.135 

0.059 

4,0,  1 

10.70 

0.13 

0.06 

4.41F-4* 

1.022 

46.3 

0.185 

1.3£-  3* 

3.  3,  1 

12.88 

0.24 

0.13 

l.6£-2 

1.061 

17.3 

0.223 

0.093 

4.4,  1 

17.48 

0.78 

0.49 

1.9£  —  1 

1.035 

29.9 

0.300 

1.112 

5,  5,  1 

19.59 

0.79 

0.48 

1.5  E-  1 

1.022 

46.5 

0.33S 

0.779 

6,6,2 

12.11 

0.13 

0.06 

9.2 E  -  6 

1.030 

34.7 

0.210 

3.6£  —  5 

7, 7.2 

14.95 

0.75 

0.50 

l.tE-2 

1.022 

45.6 

0.258 

0.324 

*£  —  I  and  E  —  3  are  factors  of  10" 1  and  10  respectively,  etc. 


tion  of  the  wave  cany  an  energy  EM.  The  quantity 
(E^/Em  )in  will  be  denoted  by  the  symbols  7f  and  Ru 
where  i  and  j  denote  the  type  of  the  incident  and  outgoing 
wave,  respectively.  The  following  subscripts  are  used  for  i 
and  j:  W  for  an  acoustic  wave  in  water,  L  for  a  longitudinal 
wave  in  the  solid  and  S  for  a  shear  wave  in  the  solid  with  the 
particle  motion  lying  in  the  plane  formed  by  the  ingoing  and 
outgoing  rays.  The  Tv  and  RtJ  defined  in  this  way  are  trans¬ 
mission  and  reflection  coefficients  which  take  on  real  values 
between  0  and  1;  phase  changes  will  be  accounted  for  by  a 
separate  phase  factor.  The  T0  and  Rtj  apply  to  an  unbounded 
plane  wave  so  that  complications  related  to  the  beam  shifts 
near  critical  angles  (e.g.,  the  Schoch  effect  and  Rayleigh 
waves)  are  not  explicitly  described  in  this  approximation. 

Expressions  for  the  TtJ  and  R,t  have  been  derived  by 
using  boundary  conditions  which  approximate  the  water  as 
an  inviscid  liquid.  These  expressions  are  well  known  in  seis¬ 
mology  (see,  e.g.,  Ref.  17,  Eq.  7  to  17  of  Ref.  18,  and  the 
review  in  Ref.  19).  The  fortran  computer  program  listed  in 
Ref.  19  facilitates  their  computation  with  minor  modifica¬ 
tions  which  make  it  possible  to  set  the  shear  wave  speed  of 
one  media  (the  water)  to  zero.  Results  from  our  version  of  the 
program  were  checked  and  found  to  agree  with  standard 
results’® 20  over  the  full  range  of  6,  the  ray  angle  relative  to 
the  normal  in  the  water.  (Reference  20  gives  Aww  ,RLL,  and 
Ass  for  an  aluminum-water  interface.)  Energy  conserva¬ 
tion,1®  symmetry  arguments,  and  numerical  computations 
verify  the  following  relationships 


TVs  —  Tsw,  rWL  —  Tlw>  7? lx  —  A SL ,  (5) 

R  LL  +  R  LS  +  T"lw  =  1, 

A  Is  +  R  SL  +  T  xw  —  1 , 

R  WW  +  T2WS  -}-  T  Vi,  =  1, 

where  the  ray  angles  6,  vs,  and  vL  are  related  by  the  general¬ 
ized  Snell’s  law.  There  are  critical  angles  of  incidence  0CJ 
such  that  7V,  =  0  for  0>0CJ,  where 

sin  0cl=Ml ,  sin  6cS=Ms.  (6) 

A  water-fused  silica  interface  at  22  'C  has  0cL  =  14.54'  and 
0cS  =  23.31°.  For  this  case  ,  RLL ,  ASL ,  and  relevant  pro¬ 
ducts  with  Twl  and  Tws  are  shown  in  Fig.  5.  Zeros  are 
evident  in  Rl  L  near  both  9.8'  and  O'L  and  in  Ass  near  both 
11.2'  and  0cL .  These  features  in  A,  ,  and  Ass  are  discussed  in 
detail  in  Ref.  18  (see,  e.g..  Figs.  1  and  7). 

To  model  the  virtual  source  strength,  the  product  of 
(£■„„,/£■,„) 1/2  at  each  vertex  is  required.  This  may  be  written 

B  —  TVah  Tii,  +  d**'  JJ  AWWJl|,  (7) 

where:  i(k  )  and  j{k )  are  the  incident  and  outgoing  wave  types 
at  the  k  th  vertex;  k  =  1  corresponds  to  the  first  refraction  of 
the  incident  ray;  and  i{k )  =j{k  —  1 ).  We  have  computed  B 
for  several  mixed  and  unmixed  chords.  Results  for  unmixed 
chord  cases  are  summarized  in  l  able  I.  Of  these  cases,  the 
largest  values  of  B  are  for  those  with  both  m  =  n  and 

en.m.l>0cL- 


TABLE  II.  Reflection  »nd  focal  parameters  of  chord  sequences  for  selected  mixed-chord  glories.  The  sphere  is  fused  silica  in  fresh  water  at  22  "C. 


(a,  m,  /( 

m.  l(4*8) 

b/a 

B 

Chord  sequence 

3, 1. 1 

8.88 

0.154 

2.1£  -  2 

LLS,  SLL 

8.88 

0.154 

1.8£  —  l 

LSL 

3,2.  1 

10.41 

0.181 

2.  If  —  2 

SSL,  LSS 

10.41 

0.181 

l.5£  -  1 

SLS 

4.2, 1 

12.84 

0.222 

l.6£  —  2 

LLSS,  SSLL 

12.84 

0.222 

9.9£  -  2 

LSLS,  SLSL 

12.84 

0.222 

1.9£-  2 

LSSL 

12.84 

0.222 

8.2£  -  2 

SLLS 

5, 3. 1 

14.45 

0.250 

2.2£  -  2 

LSLSS,  LSSLS,  SLSSL,  SSLSL 

14.45 

0.250 

2.8£  —  3 

LLSSS,  SSSLL 

14.43 

0.250 

l.$£  —  2 

SLLSS,  SSLLS 

14.45 

0.250 

1.2£  -  1 

SLSLS 

14.45 

0.250 

4.0£  —  3 

LSSSL 

806 
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For  mixed-chord  glory  rays  having  /  =  1,2  and  n  ~  3, 
4,  and  5,  B  was  computed  for  all  «!/[(«  —  nt)\m\]  chord  se¬ 
quences.  Table  II  gives  results  for  some  of  the  rays.  Evident¬ 
ly  B  depends  on  the  sequence. :i  Sequences  with  mode  con¬ 
version  at  each  reflection  tend  to  have  larger  B  than 
sequences  for  which  the  L  chords  are  grouped.  For  example, 
the  largest  B  values  of  the  (4, 2, 1 )  glory  are  for  the  LSLS  and 
SLSL  sequences;  for  the  (S,  3,  1)  glory  it  is  the  SLSLS  se¬ 
quence.  These  enhancements  of  B  are  due  to  the  relative 
magnitudes  of  RSL ,  RLL,  and  Ris  for  the  relevant  ray  angle 
in  the  water  B„m,  ■  The  equality  of  the  B  for  certain  of  the 
sequences  for  a  given  ( n ,  m,  l )  can  be  understood  via  Eq.  (5). 


III.  AMPLITUDE  AND  PHASE  AT  THE  EXIT  PLANE 

Let  Pi  exp|  —  itot )  denote  the  incident  pressure  in  the 
dashed  plane  through  C '  in  Fig.  2  where  a>  and  p,  are  the 
frequency  and  amplitude  of  the  wave.  To  an  observer  in  the 
water,  each  resulting  virtual  source  radiates  a  toroidal 
wave.9  The  amplitude  in  the  exit  plane  of  this  outgoing  wave 
is  given  by  ray  optics  to  be1 ' 


p,B 

P(s)  =  -TTT-exp 

q 


‘(vn.m.i  +  P  -  V'W  +  k{S2at> '  ) 


(81 


where  s  is  the  radius  in  the  exit  plane  and  the  parameter  p 
and  function  <F (s)  depend  on  n,  m, /as  doq,  B,  b.  and  a.  The 
phase  shift  p  accounts  for  the  phase  advance  of  rr/2  associat¬ 
ed  with  each  crossing  of  a  focal  curve  prior  to  reaching  the 
exit  plane."  Crossings  occur  at  points  L ,  and  L,  in  Fig.  2. 
The  total  p  becomes  —  ir/2  times  the  number  of  crossings 
for  a  given  | n ,  m,  I )  and  ray  sequence.  The  modulus  factor  B 
and  the  spreading  factor  q  are  each  evaluated  for  rays  with 
0=  Eq.  (8)  neglects  their  dependence  on  sin  anticipa¬ 
tion  of  our  use  of  the  stationary  phase  approximation. 

Phfse  changes  due  to  transmission  and  reflection  at 
each  vertex  are  combined  to  give  <P.  When  0n  m  l  <  Bcl ,  is  0 
or  ? t  and  is  independent  of  s.  For  unmixed  shear  glories  hav¬ 
ing  Qn.m.i  >  0cL ,  the  phase  shift  at  each  vertex  differs  from  0 
or  it;  it  depends  on  6  in  the  corresponding  plane  interface 
problem.19  Consequently  IF  depends  on  s  according  to  the 
original  angle  of  incidence  6  for  that  ray  which  crosses  the 
exit  plane  at  s.  This  dependence  is  approximated  in  Appen¬ 
dix  A  with  the  following  results  when  ka>  1 .  (i)  IF  (s|  may  be 
replaced  by  IF  (b )  in  Eq.  (8)  with  b  and  a  replaced  by  an  effec¬ 
tive  focal  radius  b  and  distance  a,.  (ii)  The  masnitudes  of 
[b  —  b,)/a  and  (a  —  a,)/a  are  roughly  (Afl)-  \  (iii)  The  cor¬ 
rections  to  b  and  a  are  so  small  in  our  experiment  that  they 
will  be  omitted.  In  the  remainder  of  this  paper  (with  the 
exception  of  Appendix  A)  IF  will  denote  IF  (6 ). 

For  the  mixed-chord  glories,  Eq.  (8)  describes  the  am¬ 
plitude  due  to  each  sequence.  The  total  amplitude  for  a  given 
(n,  m,  l )  requires  that/)'  be  summer  overall  chord  sequences. 
As  noted  in  Sec.  I,  the  a  and  q  have  not  been  modeled  for 
mixed-chord  glories. 

The  definition  of  B  as  the  product  of  (EM  /E,„  )'n  at 
each  interface  is  such  that  p'  is  specified  only  for  the  same 
medium  (water)  in  which  p,  is  given.  Modeling  of  stresses 
within  the  sphere  requires  the  use  of  material  dependent  con¬ 
version  factors  and  it  need  not  be  described  here. 


IV.  FARFIELD  SCATTERING 

In  this  section,  the  amplitude  in  the  region  ryka2  is 
considered  where  r  is  the  distance  from  C'  in  Fig.  2  to  the 
observation  point  (not  shown)  to  be  denoted  as  Q.  Let  the  left 
extension  of  the  CC '  axis  make  an  angle  y  with  the  line  C  Q. 
The  Fraunhofer  approximation  expresses  the  farfield  dif¬ 
fraction  due  to p'(s )  as  a  Hankel  transform  of  p'(s).  The  trans¬ 
form  may  be  evaluated  using  stationary-phase  and  angular- 
spectrum  techniques.13  The  farfield  pressure  due  to  a  given 
(n,m,l)  glory  ray  and  chord  sequence  will  be  denoted  by  p„ . 
The  results  in  Sec.  IIIB  of  Ref.  13  become  for  the  elastic 


sphere  case 

pn(r,r)  =  \planr\e‘'k'*"Gn,  (9) 

G„{ka,un)  =  ( ka),,2#BJlt(u„ ) 

Xexp[/(7?„m,,  +p  —  (F  — Jjt)],  (10) 
un  =  kb  sin  y,  (11) 

p  =  —  ka(  1  —  cos  y\,  (12) 

t  —  2b(2ir\a\/q)'  2a  '  :.  (13) 


w  here  for  notational  convenience  the  grouping  of  factors  in 
(9|  and  1 10)  differs  from  that  in  I F28)  and  (F29).  The  signifi¬ 
cance  of  G„ .  u„,<p.  and  /  are  as  follows:  G„  expttip)  is  the 
farfield  form  function  which  describes  the  glory  amplitude 
relative  to  the  geometric  reflection  from  a  fixed-rigid  sphere 
of  the  same  size;  u„ .  the  argument  of  the  zero-order  Bessel 
function  y„  corresponds  to  the  stationary  phase  point  of  the 
diffraction  integral  (F24)  w  hen  y  is  small;  the  phase  shift  <p  is 
due  to  a  shift  of  the  stationary-phase  point  from  its  small-y 
value:  and  t  is  a  dimensionless  geometric  parameter  which 
includes  the  effects  of  q.  b.  and  a  on  the  scatter  amplitude. 
The  use  of  the  stationary-phase  approximation  in  the  diffrac¬ 
tion  integral  motivated  our  approximation  of  B  and  q  in  Eq. 
(81  as  constants  evaluated  at  the  stationary  value  of  s  which  is 

5  -  b. 

In  addition  to  the  glory  contnbutions,  axial  rays  also 
contribute  to  the  farfield  scattering.  The  amplitudes  have 
the  following  form  according  to  ray  acoustics13 
Pn  =  Pi  (o/2r)exp(ikr  +  /£„  \f„ ,  where  n  =  0,2,4...  is  the  num¬ 
ber  of  internal  chords  and  the /„  are  the  axial  form-function 
moduli.  The  strongest  amplitude  is  due  to  the  first  axial  re¬ 
flection  (n  =  0)  which  has  /„  =  Aww  and  £0  given  by  Eq. 
(F32).  Here  /?ww  is  the  reflection  coefficient  evaluated  with 

6  =  r/2.  For  small  values  of  y,  axial  ravs  due  to  shear  chords 
are  weak  and  only  the  longitudinal-chord /„,  n  —  2, 4,...,  are 
significant.  Geometric  approximations  for  these/,  and  an 
approximation  for  follow  from  results  in  Sec.  IIC  and 
Appendix  C  of  Ref.  13. 

In  the  geometric  approximation  for  the  axial  form  func¬ 
tion  moduli,  the  /„  are  independent  of  ka.  Consequently,  if 
attenuation  is  negligible,  as  assumed  in  the  derivation  of  ( 1 0), 
glory  amplitudes  are  enhanced  by  a  factor  proportional  to 
(ka) 1 ' 2  relative  to  those  due  to  axial  rays.  When  y  =  O',  R  ww 
is  given  by  the  usual  result  for  waves  at  normal  incidence, 
Rww  =  (ZL  —  ZW)/(ZL  +  Zw),  where  the  impedances 
ZL  =  p,cL,  Zw  =  p<fw,  and p,  and p0  are  the  densities  of 
fused  silica  and  water,  respectively.  At  22  *C,  the  resulting/, 
is  0.796.  Comparison  with  the  |G„  (457,0) |  listed  in  Table  I 
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shows  that  the  glory  amplitudes  exceed  that  due  to  the  stron¬ 
gest  axial  ray. 


V.  NEARFIELD  SCATTERING 


In  this  section  we  extend  our  approximation  so  as  to 
include  points  with  r>a  which  did  not  satisfy  the  previous 
condition  r>/fca2.  It  is  again  assumed  that  ka>  1 .  The  coordi¬ 
nate  system  to  be  used  is  shown  in  Fig.  6.  The  location  of  the 
origin  C„  and  the  longitudinal  coordinate  z„  of  the  observa¬ 
tion  point  Q  depends  on  the  ray  parameters  (n,  m,  I );  C„  is  the 
center  of  the  focal  circle  of  the  (n,  m,  l )  ray.  Figure  2  shows 
one  of  these  centers.  With  z  defined  as  in  Fig.  1,  z„  =  a  4-  z. 
Scalar  diffraction  theory22,23  gives  the  following  diffraction 
integral  for  the  pressure  p„  at  Q  when  r'„  >A. 


pAxjjJ = 


Dneik'S\s-b) 


dx'  dy', 


(14) 


where  s  =  (jc'2  +  y'2)1  /2,  r'  is  the  distance  from  a  point  on  the 
focal-circle’s  plane  to  Q,  and  D„  describes  the  strength  of 
virtual  source.  Three  approximations  simplify  the  evalua¬ 
tion  of  (14).  In  the  exponent  take 


r;^z„  +  ((x  -x')2  +  l y -/)2]/2z„. 

This  is  the  Fresnel  approximation 23  which  is  applicable  pro¬ 
vided22 

^-•[(x-xf  +  ly-/)2)^.,, 

where  x'2  +  y'2  =  b 2.  Let  h  denote  the  radial  coordinate  of 
Q '  in  Fig.  6.  The  condition  becomes  zi„>X  ~\h  +  b)*  which 
for  on-axis  locations  of  the  hydrophone  reduces  to 
z„ >b  (A  /A  ),/3.  The  following  approximations  are  less  re¬ 
strictive  and  require  that  z„  >b  and  z„  >h.  Replace  r'n  in  the 
denominator  by  z„ ,  a  constant.  Let  the  focal  circle  be  uni¬ 
form  in  strength24  so  that  D„  does  not  depend  on  the  location 


y' 


FIG.  6.  Coordinate  system  which  facibutet  the  descnpoon  of  nearfleld 
scattering.  The  x'jf  plane  contains  the  focal  circle  appropriate  to  the  glory 
wave  under  consideration  and  C.  is  the  center  of  the  circle. 


of  Q  The  integral  can  now  be  simplified  by  introducing  the 
following  polar  coordinates:  x'  =  s  cos  tf>,  y'  =  s  sin 
x  =  h  cos  {■,  and  y  =  h  sin  £.  Equation  (14)  becomes 


Pn  = 


DJ 


.2w 


iiz „ 


Jflv  /»oo 

.  d*L sis 


xe,p(-  )»_»),  us. 


where  H  =  h  2/lzl„ .  This  expression  reduces  to 
Adw.)  =  (kbDJizMkbh/z^  +  (16) 


where  K  —  b  2/2zi„  and  h/z„  =  tan  y„ .  For  points  Q  not  far 
from  the  axis,  z„  ( 1  +H  (~r„ ,  the  distance  from  C„  to  Q.  To 
evaluate  D„  we  require  that  as  K~* 0,  p„(0, 

\-*pn  ( y  =  0,  r  —  zn  —  a)  from  Eq.  (9).  This  evaluation 
when  combined  with  ( 1 0)  yields  the  following  approximation 
for  the  pressure  in  the  near-  and  farfields 

Pn  )  =  {p,a/2zn  )exp  [  ik  (r„  -  a  +  z,K )  ] 

xGn\ka,kbh/zn),  (17) 


where  G„  is  given  by  ( 10)  evaluated  with  uH  =  kbh  /z„ .  It  can 
be  shown  that  when  y  is  not  large  and  tykb2  that 
k  (r„  —  a  4-  znK  )-*{kr  +  <p ).  Consequently  (17)  reproduces 
the  previous  farfield  result  both  on  and  off  of  the  axis. 

To  approximate  the  pressure  due  to  the  axial  reflection 
in  the  near  field,  we  use  the  geometric  result13  that  for  h  / 
z<  1 ,  the  reflection  appears  to  come  from  a  virtual  pointlike 
source  located  a  distance  a/2  behind  C This  source  is  at  A0 
in  Fig.  2  and  Ref.  13,  Fig.  4.  The  reflected  pressure  at  Q 
becomes 


F0  =  (p,a/2r0)Rwwe'i'f»~2  (18) 

where  r„  is  the  distance  from  A„  to  Q  and  R  ww  is  evaluated  at 
the  geometric  angle  of  incidence.  From  Weston’s  asymptotic 
results  for  electromagnetic  nearfield  backscattering,”  it  is 
expected  that  the  leading  correction  to  ( 1 8)  will  be  of  magni¬ 
tude  {ka)~ 1  (or  smaller)  relative  to  unity  provided  r0  is  some¬ 
what  larger  than  a/2. 


VI.  TRANSIENT  PULSES:  DELAY  AND  DISTORTION 


In  this  section  the  pressure  of  the  incident  plane  wave  is 
changed  from  p,ex p(  —  ioit )  to  a  transient  function  of  time 
p,(t )  where  the  reference  plane  is  again  the  dashed  plane  in 
Fig.  2.  To  simplify  our  description  of  the  nearfield  scattering, 
Q  is  now  restricted  to  lie  on  the  backward  axis  so  that 
x  =  y  =  0.  After  the  incident  wave  crosses  the  dashed  plane, 
the  propagation  time  delay  r  required  to  reach  Q  is  deter¬ 
mined  by  the  combined  k -dependent  phase  terms  in  (10)  and 
(17): 

r  =  r0  +  (19) 

cwTn.m.,  =  2a  [  1  -  cos  eH  m  ,  +  mMs  cos  vs 

+  (n  -m)ML  cosvt]  -f^^-ha)]-'.  (20) 
r0  —  z/cw  is  the  propagation  delay  for  the  first  axial  reflec¬ 
tion  to  reach  Q.  For  the  [n,  m,  I ) th  glory  wave,  r,  _  ,  is  the 
delay  relative  to  that  of  the  axial  reflection.  The  pressure  of 
the  first  axial  reflection  is 

pS )  =  affww(2z  +  a)~  'pr(t  -  r0}. 

Unlike  the  axial  wave,  the  glory  waves  are  distorted  as 
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well  as  shifted  in  time.  To  describe  the  distortion,  it  is  con¬ 
venient  to  define 

e=n-*,  (21) 


where  p  and  V  are  phase  shifts  due  to  internal  caustics  and 
reflection-transmission  factors,  respectively  (Sec.  III).  Since 
e  <  0,  the  phase  is  advanced  and  this  ^-independent  shift  dis¬ 
torts  transients.23-26,27  Further  distortion  is  predicted  as  a 
consequence  of  the  k  1/2exp(  —  <ir/4)  factor  in  (10).  Appendix 
B  models  the  combined  distortions;  combined  with  (17)  and 
( 19)  it  yields  the  following  for  the  pressure  of  the  (n,  m,  /  )th 
glory  wave: 


/>*(')  = 


ainWB  f  d  \ 
2z„cl?  I  dt ) 


1/2 

PAt-r), 


(22) 


where  [d /dt  )l/2  is  the  half-order  derivative  operator  [Eq. 
(BS)],  p,  is  given  by  (B3)  evaluated  at  (t  —  r ).  In  the  evalua¬ 
tion  of  (B3 ),pH(t  —  r)  is  the  Hilbert  transform26  of p,{t  —  r) 
which  follows  from  (B4). 

The  phase  shift  by  e  of  p,{t  —  r)  yields  p,(t  —  r).  This 
description  of  the  distortion  is  only  approximate  since  the 
shift  e  is  only  applicable  to  spectral  components  having 
ka>  1.  Application  of  this  result  to  p,(t )  which  turn  on  ab¬ 
ruptly  at  t  —  0  predicts  precursor  in  pH  which  extends  un- 
physically27  to  t  =  ~  oo .  Due  to  the  asymptotic  limitations 
of  (17),  the  resulting  temporal  expressions  (B3)  and  (22)  in¬ 
correctly  describe  the  early  time  behavior  of  the  precursor. 
In  the  experiment  to  be  described  the  incident  pressure  is 
approximated  by 


p,[t  I  =  Pi  sin(<ur ),  0<t  <2irN /cv,  (23) 

with  p,(t )  =  0  elsewhere;  M  is  an  integer  number  of  cycles, 
each  of  amplitude  p,.  Cron  and  Nuttall26  have  graphed  the 
resulting  pf(t )  for  /V  =  1  and  2  and  various  e  (their  Figs.  2 
and  4).  Their  analysis  indicates  that  for  integer  N>  2,  the 
peak-to-peak  amplitude  of  pt  in  the  central  cycles  of  the  tone 
burst  is  close  to  the  undistorted  value,  2p,.  Furthermore  for 
these  A,  both  the  after-signal  and  the  precursor  are  weak. 

To  estimate  the  effect  of  the  half-order  derivative  on  the 
peak-to-peak  amplitude  of  the  scattering  we  have  considered 
( d  /dt)' 12  p ,  forthep,  given  by  (23).  When  A>4,  the  peak-to- 
peak  amplitude  of  the  central  cycles  of  the  resulting  function 
is  within  0.2%  of  the  steady-state  value  2|<u|  U2p,  which  fol¬ 
lows  from  (B6).  This  steady-state  approximation  applies  to 
the  scattered  waveform  when  sin  e  =  0;  however,  it  is  to  be 
expected  that  even  when  sin  e^G,  it  approximates  ihe  peak- 
to-peak  {d /dt  )'r2p,  due  to  the  aforementioned  smallness  of 
the  effect  of  e  on  the  central  peak-to-peak  amplitude.  The 
effect  of  the  id  /dt  )l/2  operator  on  certain  other  functions  has 
been  tabulated.28  2g 

In  addition  to  the  distortions  due  to  e  and  the  (d  /dt  )1/2 
operator,  broadbanded  signals  in  water  are  distorted  due  tc 
the  preferential  attenuation  of  the  high  <y  components.  This 
distortion  is  estimated  to  be  negligible  in  the  experiment  to 
be  described. 


VII.  EXPERIMENTAL  CONSIDERATIONS  AND 
HYDROPHONE  CALIBRATION 

This  section  supplements  the  experiment's  description 
given  in  the  Introduction.  We  note  considerations  which  fa¬ 


cilitate  a  test  of  the  model  developed  in  Secs.  V  and  VI.  Fig¬ 
ure  7  shows  the  scattering  chamber  (an  8X2X2$  ft  water 
aquarium),  the  transducers,  and  the  sphere.  The  sphere  was 
supported  by  three  plastic-coated  aluminum  rods.  (The  usu¬ 
al  practice13  of  hanging  the  target  from  a  wire  was  not  used 
due  to  risks  of  damage  to  the  expensive  sphere.)  The  rods 
were  positioned  so  as  to  avoid  paths  of  glory  rays.  Rod  ech¬ 
oes  were  identified  and  eliminated  using  strips  of  anechoic 
rubber  (not  shown  in  Fig.  7).  Echoes  from  the  tank's  rear  wall 
were  similarly  eliminated.  The  water  was  filtered  and  deion¬ 
ized.  The  x  value  of  the  hydrophone  probe  (Fig.  1 )  was  con¬ 
trolled  by  the  stepper  motor  which  is  visible  in  Fig.  7.  Rails, 
which  extended  the  length  of  the  tank,  supported  the  driven 
transducer  and  the  hydrophone  positioner.  The  rails  facili¬ 
tated  the  location  of  the  backward  axis  (the  truex  =  O.y  =  0, 
location)  to  within  1  mm.  Distances  zp  and  z  were  deter¬ 
mined  sonically  by  measuring  tone-burst  propagation  delays 
with  the  aid  of  a  digital  delayed-pulse  generator. 

The  hydrophone’s  piezoelectric  element  was  sufficient¬ 
ly  small  (0.6  mm  diam)  that  its  directivity  was  insignificant 
for  the  range  of  x  encountered  in  the  experiment.  The  bent 
hypodermic  needle  which  supported  the  hydrophone  was 
sufficiently  small  in  diameter  (1.3  mm)  that  its  perturbing  of 
the  incident  wave  was  insignificant  in  the  measurements  of 
the  glory  scattering. 

The  source  transducer  (Panametrics  model  V309,  ac¬ 
tive  element  diameter  of  d  =  12.7  mm)  was  excited  at  fre¬ 
quencies  below  its  S-MHz  resonance.  It  was  driven  by  a  low 
impedance  source  consisting  of  two  matched  linear  40-W 
power  amplifiers  (Amplifier  Research  model  40A12)  operat¬ 
ed  with  their  outputs  in  parallel.  The  voltage  V (/ )  across  the 
transducer  was  closely  approximated  by  the  function 
V  —  V  sin(oit )  for  0  <  f  <  2irN  /a>  and  P  =  0  all  other  t,  where 
N  =  4.  This  excitation  was  repeated  at  ~  10-ms  intervals; 
this  interval  was  sufficiently  long  for  the  decay  of  spurious 
reflections.  The  tranducer’s  radiation  pattern  with  to/ 
2 n  =  2. 1  MHz  was  measured  near  the  z  =  0  plane  as  a  func¬ 
tion  of  the  impact  parameter  b.  (This  was  facilitated  by  re¬ 
versing  the  direction  of  the  hydrophone  from  that  shown  in 
Figs.  1  and  7.)  Measurements  of  the  incident  wave’s  ampli¬ 
tude)?,  as  a  function  of  b  (over  the  range  of  relevant  b )  agreed 


FIG  7.  Photograph  of  the  scattering  facility.  The  source  transducer  is  on 
the  left,  the  bent  hydrophone  needle  is  in  the  center,  and  the  sphere  is  on  the 
right.  Adjacent  to  the  base  of  the  aquarium  is  a  rule  marked  in  centimeters. 


J.  Accost.  Soc  Am.,  Vol.  74,  No.  2,  August  1983 


Marston  eta/. .  Back  sea  tiering  from  sphere 


611 


611 


approximately  with  the  elementary  diffraction  model  of  the 
central  radiation  pattern  of  a  piston:  p,(b ) 
~p,(6  =  0)(  1  —  b  2v)  where10  v  =  \kd  /2 zp  )2/8.  Glory  rays  in 
our  model13  have  b<bm,k  =  4/Aa~2.0  cm.  The  resulting 
correction  to  Eq.  (17)  is  small  yet  significant:  b  ^,xi'~0.05;  it 
is  incorporated  into  Eq.  (24). 

Figure  8  shows  oscilloscope  records  of  the  first  axial 
echo  and  several  glory  echoes.  Figure  8(c)  snows  the  axial 
echo,  and  after  delays  r„  m/  >  40  /is,  several  distinct  glory 
echoes.  Inspection  of  the  time  interval  between  the  axial  and 
glory  echoes  indicates  that  the  noise  amplitude  is  ~  10%  of 
the  typical  echo’s  amplitude.  Tests  indicate  that  the  princi¬ 
pal  noise  source  was  the  wideband  preamplifier  which  was 
matched  to  the  hydrophone.  The  preamplifier's  self-noise 
was  equivalent  to  an  input  ~2  p V  rms;  however,  the  hydro¬ 
phone  output  voltage  was  small  due  to  the  combined  require¬ 
ments  of  high  spatial  resolution  (small  size  of  piezoelectric 
element)  and  good  temporal  resolution  (heavily  damped  ele¬ 
ment  for  wide  bandwidth).  Not  having  signal  averaging  ca¬ 
pabilities  limits  the  extent  of  these  experiments  to  measure¬ 
ments  of  echo  times  and  (as  a  function  of  x)  amplitudes  but 
not  their  distortions. 

Our  procedure  of  testing  Eq.  ( 1 7)  is  to  first  measure  the 
peak-to-peak  preamplifier  output  voltage  for  the  central  cy¬ 
cles  of  the  axial  echo.  This  voltage  is  measured  when  x  =  0;  it 


(c) 

-•  4  sftb 

10  ps 

FIG.  8.  Photographs  of  oscilloscope  traces  of  the  preamplifier  output  for 
the  hydrophone  on  axis  and  r~36  cm.  (a)  Shows  the  first  axial  reflection,  (b) 
Begins  67  /is  after  the  axial  reflection  and  shorn  several  glory  echoes,  (c) 
Begins  a  few  /is  before  the  reflection  and  it  shows  both  the  reflection  and  the 
principal  glory  echoes.  The  strongest  echo  in  (b)  and  (c)  is  the  (4, 4.  I |  glory 
echo. 


is  denoted  as  xP„(x  =  0)  where  k  is  an  undetermined  constant 
and  Pu  denotes  the  pressure  amplitude.  Then,  for  a  specific 
glory  echo  (which  is  identified  from  its  TnnJ),  the  peak-to- 
peak  voltage  kP(x)  is  measured  for  the  central  cycles  of  the 
echo.  The  measurements  are  made  from  the  centers  of  the 
noise  broadened  oscilloscope  traces  to  facilitate  a  partial 
cancellation  of  noise.  The  result  of  Sec.  VI  is  that  the  steady- 
state  model  is  applicable  to  the  peak-to-peak  central  portion 
of  the  glory  echo’s  amplitude.  Consequently,  the  model  pre¬ 
diction  for  the  voltage  (and  amplitude)  ratio  is 

P(x)  _  (*  +  W\Gn{ka,kbx/zn)\{\  -  b2v) 

P '„(x  =  0) 

where  the  aforementioned  source  factor  (1  —  b  2o)  is  included 
and  Rw  is  evaluated  at  normal  incidence. 

A  technique  for  calibrating  hydrophones  which  refer¬ 
ences  elastic  response  to  the  amplitude  of  the  axial  reflection 
is  described  by  Dragonette  et  a/.11  Our  technique  uses  geo¬ 
metric  and  reflection  factors  whch  differ  from  those  of  Ref. 
31  in  ways  evident  from  Eqs.  (18)  and  (24).  (i)  The  spreading 
factor  in  the  backscattered  axial  echo’s  amplitude  is 
(z  +  [a)“ ',  not  (z  +  a)  - 1  as  it  would  be  if  the  wave  appeared 
to  spread  from  the  sphere’s  center,  (ii)  The  reflected  ampli¬ 
tude  is  less  than  that  for  reflection  from  a  rigid  sphere  due  to 
the  coefficient  Rww .  That  the  amplitude  of  the  specular  re¬ 
flections  is  altered  by  a  target’s  impedance  has  been  con¬ 
firmed  previously  with  exact  calculations  of  the  backscatter- 
ing  of  tone  bursts,  e.g.,  those  of  Davis  et  al? 2  for  various 
rubber  cylinders  with  ka  —  12. 

Some  comments  on  attenuation  due  to  bulk  absorption 
of  the  2. 1  MHz  bursts  are  in  order.  The  bulk  absorptions  of 
shear  and  longitudinal  waves  within  the  fused  silica  are  cal¬ 
culated  to  be  negligible  from  published  attenuation  data 
even  w  hen  the  internal  path  length  was  long,  as  is  the  case  for 
the  (7,  7.  2)  path  shown  in  Fig.  3.  [The  reasons  for  choosing 
fused  silica  were  (il  low  attenuation  even  at  high  frequencies, 
(ii)  isotropy  of  elastic  properties,  and  (iii)  dimensional  stabil¬ 
ity  .]  Absorption  losses  within  the  water  over  the  propagation 
distances  zr  and  z  are  not  negligible  if  absolute  amplitude  is 
the  quantity  of  interest.  The  experiments,  however,  yielded 
the  ratio  of  axial  and  glory  amplitudes.  The  path  lengths  in 
water  of  the  axial  and  glory  echoes  differed  by  <  1  cm.  Con¬ 
sequently  these  losses  do  not  significantly  alter  [Eq.  (24)],  the 
model  prediction. 

The  water  temperature  T  varied  with  the  room  tem¬ 
perature.  Associated  with  periods  of  rapid  change  in  T,  the 
glory  echo  amplitudes  sometimes  contained  irregularities  in 
their  dependence  on  x.  These  irregularities  were  apparently 
due  to  the  influence  of  gradients  in  T  on  the  acoustic  propa¬ 
gation.  Data  collection  was  limited  to  periods  for  which  the 
irregularities  were  absent.  During  theses  periods,  measure¬ 
ments  of  T\y )  (via  a  vertical  array  of  thermalcouples)  gave 
vertically  averaged  gradients  (dT  /dy)  having  moduli 
<1  *C/m  though  local  gradients  were  as  large  as  1  *C/m. 

During  the  series  of  experiments,  the  mean  water  tem¬ 
perature  (T)  ranged  from  22*-25  *C.  Details  of  the  r,m,  and 
P [x)/P0  data  were  only  weakly  dependent  on  <T>;  however, 
comparisons  with  the  theory  are  always  made  for  material 
parameters  appropriate  to  the  T  of  the  measurement.  The 
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fresh  water  data  of  Ref.  33  for  were  used;  e.g.,  at  22  *C,  cw 
=  1 .488  km/s.  Use  was  made  of  Fraser’s  measurements34  of 
cs  and  cL  for  synthetic  fused  silicas  having  the  normal  OH 
ion  concentration  of  103  ppm  by  weight.  The  measurements 
were  adjusted  for  T  of  22*  to  25  *C  by  using  the  temperature 
coefficients35:  dcs/dT~0.2S  m/s  *C  and  dcL/dT~ 0.68  m/ 
s  *C.  At  22  'C,  this  procedure  gives  cs  =  3.761  km/s  and  cL 
=  5.928  km/s.  The  density  of  synthetic  fused  silica  is  2.201 
g/cm3;  at  22  *C,  the  density  of  water  is  0.9978  g/cm3. 

VIII.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

A.  On-axis  acho  amplitudes  and  dalaya 

We  now  compare  measurements  taken  from  oscillo¬ 
scope  traces,  similar  to  those  of  Fig.  8,  with  the  model  pre¬ 
dictions  of  Eqs.  (20)  and  (24).  The  hydrophone  was  placed 
within  1  mm  of  the  true  backward  axis.  (The  proximity  of 
this  axis  is  evident  from  the  symmetry  of  off-axis  amplitude 
data  described  later.)  Time  delays  were  measured  from  the 
central  cycle  of  the  axial  reflection  to  the  corresponding  cy¬ 
cle  in  the  glory  echo.  The  latter  time  could  only  be  deter¬ 
mined  with  an  uncertainty  which  slightly  exceeded  ±J  cy¬ 
cle  ( ±  0.24  (is).  The  sphere’s  radius  was  chosen  to  be 
sufficiently  large  that  most  principal  glory  echoes  would  not 
overlap  as  illustrated  by  Fig.  8.  Representative  delay  and 
amplitude  measurements  for  distinguishable  echoes  are  giv¬ 
en  in  Table  III  together  with  predictions.  The  echoes  are 
ordered  according  to  their  delay. 

The  difference  between  the  observed  and  modeled  delay 
in  each  case  was  <0.5  (is.  When  making  this  comparison  it 
should  be  remembered  that  Eq.  (20)  does  not  include  the 
effects  of  the  ^-independent  phase  shift  e  on  the  temporal 
displacement  of  the  cycles  near  the  center  of  each  echo,26  nor 
does  it  include  displacements  due  to  the (d/dt)'12  operator  in 
Eq.  (22).  We  estimate  the  magnitude  of  these  apparent  time 
shifts  to  be  S  0.25  (is.  The  agreement  of  the  model  with  these 
measurements  facilitates  the  identification  of  these  echoes. 

One  prominant  echo  appeared  to  be  the  superposition 
of  the  (5,  3,  1)  and  (4,  3,  1)  mixed  chord  echoes;  both  are 
predicted  to  have  rnm  ,~70.0  (is.  The  resulting  superposi¬ 


tion  is  the  first  echo  in  Fig.  8(b).  It  is  followed  by  large  distin¬ 
guishable  echoes  in  Fig.  8(b)  which  are  the  (4, 4, 1)  and  (5, 5, 
1 )  echoes.  These  are  followed  by  a  broadened  echo  which  is  a 
superposition  of  echoes  for  rays  having  n  —  m>  5  and  /  =  1 . 

Other  than  the  axial  reflection,  Fig.  8(a),  only  one  other 
axial-ray  echo  could  be  identified.  This  was  the  (2, 0, 1)  echo 
listed  in  Table  III  which  has  a  predicted  delay 
r201  =4 a/cL.  The  signal-to-noise  ratio  was  noticeably 
smaller  for  z  =  74.4  cm  than  for  36.4  cm;  consequently  this 
echo  was  not  clearly  detected  in  the  large  z  case. 

For  the  measurements  with  z  =  36.4  cm,  the  contribu¬ 
tion  to  rnmJ  of  the  second  term  on  the  right  of  Eq.  (20)  was  as 
large  as  0.24  (is.  Evidently  this  term  is  significant  since  its 
omission  decreases  agreement  with  the  data.  This  term  has  a 
simple  physical  interpretation.  It  approximates  the  time  dif¬ 
ference  for  propagation  along  a  straight  path  in  water  from 
the  focal  circle  to  the  hydrophone  with  that  from  the  on-axis 
point  C„  to  the  hydrophone.  When  computing  this  term  in 
the  mixed  chord  cases  the  exact  expression  for  a  is  not 
known.  We  take  a  —  a  since  elementary  considerations  indi¬ 
cate  that  the  C„  lie  close  in  the  sphere’s  center.  The  values  of 
a/ a  tabulated  in  Table  I  give  evidence  of  proximity  in  the 
case  of  unmixed  chords. 

Let  e  and  denote  the  ratio  of  peak-to-peak  ampli¬ 
tudes  P{x  =  0)/Fo(x  =  0)  as  measured  and  as  predicted  by 
Eq.  (24),  respectively.  Table  III  compares  and  with 
two  values  of  z  at  the  indicated  temperatures.  For  most  of  the 
echoes,  was  measured  on  more  than  one  day  at  the  indi¬ 
cated  conditions  and  the  value  is  an  average.  With  apparent¬ 
ly  similar  water  condition  individual  measurements  varied 
by  as  much  as  ±  0.04  from  the  indicated  This  uncer¬ 
tainty  is  also  representative  of  comparisons  among  measure¬ 
ments  taken  at  other  temperatures.  The  deviations  often  ex¬ 
ceeded  +  0.04  in  measurements  having  z  =  74.4  cm  since 
the  signal  voltage  was  smaller  relative  to  the  amplifier's 
noise.  Inspection  of  Fig.  8(a)  reveals  a  weak  signal  which 
appears  ~2  (is  after  the  first  axial  reflection.  Tests  indicate 
that  this  was  due  to  the  hydrophone.  This  signal  may  intro¬ 
duce  a  nonrandom  source  of  error  in  the  measurement  of  $  e 
when  the  echo  follows  an  earlier  echo  by  $  2  (is  as  is  the  case 


TABLE  1 1 1.  Predicted  and  measured  backscattering  delays  and  amplitudes  relative  to  those  of  the  specular  echo.  The  upper  group  was  obtained  with  z  =  36.4 
cm  and  T  =  24.9  *C.  The  lower  group  had  z  =  74.4  cm  and  T  =  22.0  “C.  The  incident  wave  had  u>/2v  =  2.10  MHz. 


Echo 

In,  m.  1 1 

Delay  r„.  „ 

Eq.  1201 

Measured 

Ratios  of  peak-to-peak  amplitudes 

2.0,  1 

34.8 

34.8 

0.049 

0.04 

0.8 

3,0,  1 

44.6 

44.3 

0.069 

0.08 

1.16 

3,  2.  1 

62.0 

62.1 

0.54 

4.2,  1 

63.4 

63.5 

0.70* 

4.4.  1 

74.8 

74.7 

1.271 

1.23 

0.97 

5.  5,  I 

77.0 

77.0 

0.883 

0.90 

1.02 

7. 7,2 

147.9 

148.1 

0.373 

0.27 

0.72 

3,0,  1 

44.6 

44.5 

0.074 

3.  1.  1 

53.5 

53.5 

0.87 

4,4.  1 

74.7 

74.7 

1.036 

1.19 

0.91 

5,  5,  1 

77.0 

76.5 

0.907 

0.87 

0.96 

7,  7,2 

147.9 

148.0 

0.384 

0.35 

0.91 

*  The  (4,  2,  1 1  amplitude  may  contain  interference  from  the  15,  2,  I )  echo;  r,  2  ,  from  Eq.  (20)  is  63.8  ft*. 
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for  (5,  5, 1). 

For  unmixed  chord  glory  echoes,  Eq.  (24)  was  used  to 
compute  0?,  and  the  ratios  0?,/0P,  were  computed.  The 
theory  overestimates  the  amplitudes  of  the  glory  echoes  with 
the  uncertain  exceptions  of  (3, 0, 1 )  and  (5, 5, 1 )  echoes.  Note, 
however,  that  the  spread  of  values  for  29 ,/29 ,  is  much 
smaller  than  the  spread  of  the  29, .  Evidently  the  theory  cor¬ 
rectly  describes  the  relative  ordering  of  glory  echo  ampli¬ 
tudes.  The  largest  echo  is  the  (4,  4.  1);  its  amplitude  on  the 
backward  axis  is  almost  as  large  as  that  for  a  perfectly  re¬ 
flecting  sphere  of  the  same  size.  In  terms  of  P0  and  Aww  of 
the  fused  silica  sphere,  the  perfect  reflector  would  have  P  / 
/>0  =  (/lwwp'~1.256.  This  value  is  only  slightly  larger 
than  the  29,  of  the  (4, 4, 1)  echo. 

The  relative  values  of  0„  mJ  and  dcL  significantly  influ¬ 
ence  the  29,  for  those  echoes  with  n  =  m.  The  (4, 4, 1),  (5,  3, 
1),  and  (7, 7,  2)  glory  rays  all  have  0„  m ,  >  0cL .  For  compari¬ 
son,  the  (3,  3,  1)  echo  is  predicted  to  be  small,  having 
21 ,  =0.11.  Equation  (20)  gives  r,  j  ,  =  69.8  jis  which  is 
close  to  the  delays  for  the  (5, 3, 1 )  and  (4, 3, 1 )  echoes.  Conse¬ 
quently  the  (3,  3, 1)  echo  could  not  be  discriminated  from  the 
(5,  5,  1  )and  (4,  3,  1)  echoes.  For  aluminum  spheres  in  water, 
0,  3  I  >  6cL  and  the  (3,  3,  1 )  echo  is  predicted  to  be  the  stron¬ 
gest  unmixed-chord  glory  echo.9 

Consider  now  the  29,  for  the  distinguishable  mixed- 
chord  glory  echoes  (3,  1, 1),  (3,  2, 1),  and  (4,  2,  1).  For  each  of 
these  cases  there  are  several  chord  sequences  as  illustrated  in 
Fig.  4.  Table  II  indicates  that  the  B  factor  for  certain  of  these 
sequences  are  similar  in  magnitude  to  the  B  of  the  stronger 
glories  having  n  =  m.  It  is  consistent  with  our  model  that 
these  three  echoes  were  easy  to  detect,  though  their  param¬ 
eters  a  and  q  are  not  known. 

Table  III  includes  39,  and  39,  for  the  (2, 0, 1 )  axial  echo. 
Here  2P,  =  (z  +-  Jo)[/'2|(z2/!ww )~ 1  where  extension  of  results 
in  Appendix  C  of  Ref.  1 3  to  the  elastic  sphere  case  with  0  =  0 
give  the  following:  \f2\  —  Af/Iww(l  —  R  Ul,  M  =ML/ 
(2  -ML),z2  =  z  +  a2,  and  a2/a  =  1  +  \M.  The  resulting 
39 ,  agrees  with  the  measured  29 ,  within  experimental  error. 

B.  Off-axis  echo  amplitudes 

The  dependences  of  the  peak-to-peak  glory  amplitudes 
on  x  (the  displacement  from  the  backward  axis)  give  direct 
evidence  of  axial  focusing.  First  consider  measurements  of 
the  peak-to-peak  voltage  kP0  for  the  simply  reflected  [n  =  0) 


echo.  The  reflected  wave  is  not  axially  focused.  For  both  z  of 
36.4  and  74.4  cm,  there  was  no  measurable  dependence  of 
kP0  on  x  for  displacements  |x|  5  2  cm.  In  sharp  contrast  with 
this  behavior,  the  peak-to-peak  glory  echo  voltages  kP  were 
strongly  dependent  on  x. 

Figure  9  displays  the  measured  ratio  P  (x)/P^x  =  0)  as¬ 
sociated  with  the  (4  4,1)  glory  when  z  =  36.4  cm.  The  proce¬ 
dure  which  was  followed  was  to  record  kP  after  each  incre¬ 
ment  in  x.  Starting  with  x  =  0,  the  scan  was  completed  in 
one  direction;  x  was  reset  to  0  and  the  scan  completed  in  the 
other  direction.  The  increments  were  uniform;  in  most  scans 
they  were  250 fim.  Division  of  the  data  by  xP^x  —  0)  yielded 
the  plotted  points.  With  increasing  x,  the  echo's  delay  in¬ 
creases  due  to  the  x  dependence  of  r„  in  Eq.  (17). 

To  find  if  the  dependences  on  x  is  described  by  the 
steady-state  model,  the  following  quantity  was  graphed  to¬ 
gether  with  each  data  set: 

29,  =  29,\J0[kbx/z„)\,  (25) 

where  29 ,  is  the  measured  P  (0)/P„(0).  Thex  dependent  factor 
is  that  of  Eq.  (24).  Normalization  to  29 ,  was  used  since  in 
plotted  cases  29,  <  29,  where  29,  is  given  by  (24)  with  x  =  0. 
It  is  evident  from  Fig.  9  that  (25)  approximates  the  data  ex¬ 
cept  near  the  minima. 

Most  of  the  scans  were  taken  with  z  =  74.4  cm  where 
the  maxima  were  sufficiently  broad  that  the  finite  width  of 
the  hydrophone’s  active  element  (0.6  mm)  was  insignificant. 
Figures  10  and  1 1  give  results  for  unmixed-chord  and  mixed- 
chord  echoes,  respectively.  It  is  again  evident  that  the  echoes 
are  localized  or  “focused”  along  the  backward  axis.  Fluctu¬ 
ations  of  the  data  are  larger  than  in  Fig.  9  due  to  the  lower 
amplitude  of  the  voltages.  When  graphing  (25)  in  Fig.  1 1,  it 
was  necessary  to  assume  a  value  for  a  to  be  able  to  evaluate 
z„  =  z  +  a.  Here  we  take  a —  a  because  of  the  expected 
proximity  of  the  C„  to  the  sphere’s  center.  The  graphs  in 
Figs.  9  and  10  use  the  modeled  a  given  by  Eq.  (FlOb).  The 
values  of  29 ,  for  these  data  sets  are  similar  to  those  of  Table 
III. 

Each  data  set  is  the  result  of  a  single  scan  and  each  is 
plotted  without  smoothing.  Figures  9  and  10  are  representa¬ 
tive  of  1 1  scans  of  (4,  4,  1)  or  (5,  5,  1)  echoes  obtained  at 
various  water  temperatures.  Some  comments  on  the  symme¬ 
try  of  the  data  in  Figs.  10  and  1 1  are  merited.  To  symmetrize 
these  sets,  the  assumed  origin  was  shifted  by  0.5  mm  for  all 
sets  having  z  —  74.4  cm.  This  shift  is  consistent  with  the 


FIG.  9.  The  peak-to-peak  voltage  for 
the(4,4, 1)  echo  plotted  as  a  function 
of  x.  The  voltages  are  ratioed  with 
the  on-axis  value  of  the  first  axial 
echo.  Measurements  were  obtained 
when  the  water  temperature  was 
24.9  *C.  The  solid  curve.  Eq.  (25).  is 
the  modeled  dependence  on  x  nor¬ 
malized  so  as  to  fit  the  measurement 
at  x  =  0.  The  peak  at  x  =  0  and  the 
side  lobes  are  direct  evidence  of  axial 
focusing. 
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FIG.  10.  Data  and  models  (as  in  Fig. 
9)  for  the  (4, 4,  1)  and  (5,  S,  1)  echoes 
in  (a)  and  (b),  respectively.  Here 
z  =  74.4  cm  and  the  temperature 
was  22  *C. 


aforementioned  alignment  imprecision.  Since  Figs.  10  and 
1 1  were  plotted  with  a  common  x  =  0,  y  =  0  location,  the 
symmetry  of  each  data  set  is  experimental  confirmation  that 
the  different  echoes  are  focused  along  a  common  axis.  The 
predicted  dependences  on  displacement  do  not  depend  on 
the  azimuthal  angle  as  a  consequence  of  the  sphere's  isotropy 
and  placement  at  the  center  of  source's  radiation  pattern. 
The  isotropy  of  the  sphere  was  confirmed  by  comparing  (4, 4, 
1)  measurements,  Fig.  12,  obtained  before  and  after  rotation 
of  the  sphere  by  an  angle  ~  1 30*. 

Figures  9-1 1  show  minima  which  are  not  as  deep  as 
predicted  by  the  steady-state  theory.  Recall  the  assumption 
of  x  =  0  and y  =  0  in  the  derivation  of  Eq.  (22)  and  the  ensu¬ 
ing  discussion  of  the  applicability  of  steady-state  theory.  For 
off-axis  points,  the  Bessel  function  in  Eq.  ( 10)  contains  a  de¬ 
pendence  on  co  not  included  in  the  derivation  of  (22).  It  is 
assumed  in  Eqs.  (24)  and  (25)  that  steady-state  peak  to- peak 
amplitudes  are  applicable  to  off-axis  points  provided  the  inci¬ 


dent  wave  is  described  by  Eq.  (23)  with  N  —  4.  This  assump¬ 
tion  is  suspect  for  x  near  the  zeros  of  J0(kbx/zm )  where  the 
dependence  on  to  is  most  significant.  Careful  observations  of 
the  waveform  near  the  temporal  center  of  the  echo  indicate 
an  apparent  phase  shift  by  n  as  the  hydrophone  is  scanned 
through  each  minima.  This  is  in  agreement  with  steady-state 
theory.  Unfortunately,  the  peak-to-peak  amplitudes  are  dif¬ 
ficult  to  measure  near  these  minima  due  to  a  small  signal-to- 
noise  ratio. 


IX.  CONCLUSIONS  AND  RELEVANCE  OF  PARTIAL- 
WAVE  THEORY 

We  observe  that  the  backscattering  of  sound  from  elas¬ 
tic  spheres  manifests  a  localized  enhancement  similar  to  that 
present  in  the  optical  glories  of  drops14  ’6  and  bubbles. 15  The 
scattering  associated  with  particular  internal  rays  is  ob- 


FIG.  11.  Data  and  model  (at  in  Fig. 
9)  for  the  (3, 1, 1 )  mixed-chord  echo. 
Here  z  »  74.4  cm  and  the  tempera¬ 
ture  was  24  *C. 
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TRANSVERSE  DISPLACEMENT  *  OF  HYDROPHONE  (cm) 

FIG.  12.  Data  and  model  (as  in  Fig.  9)  but  fori  =  74.4  cm  and  a  temperature 
of  24.7 'C.  Data  was  obtained  before  and  after  rotation  of  the  sphere 
through  an  angle  ~  1 30*.  The  rotation  axis  was  inclined  to  the  backward 
axis;  consequently  the  rays  traversed  regions  of  the  sphere  which  differed 
before  and  after  the  rotation. 


served  to  be  axially  focused.  Inspection  of  Figs.  9-12  reveal 
angular  widths  Ay  of  a  few  degrees  for  the  central  peaks;  for 
example  in  10(a),  d/~2.8  cm/z~38  mr~2.2‘.  The  agree¬ 
ment  with  the  modeled  focusing  pattern  and  the  form  of  Eqs. 

( 10)  and  (11)  suggest  that  d/~u/ 1  /vb  radians  where  w~2.4 
is  the  first  root  of  J„[u)  =  0.  This  approximation  should  be 
applicable  when  the  hydrophone's  distance  z  and  sphere’s 
radius  a  are  such  that  z>a>> l.  Inspection  of  Figs.  10  and  1 1 
together  with  Tables  I  and  II  demonstrate  that  the  observed 
width  of  the  central  peak  is  smallest  for  echoes  having  the 
largest  focal  radii  b.  Inspection  of  Table  III  shows  that  the 
discrepancy  between  the  observed  and  predicted  central  am¬ 
plitude  is  approximately  5%  of  the  predicted  value.  The 
cause  of  this  discrepancy  is  not  known;  however,  experimen¬ 
tal  difficulties  due  to  preamplifier  noise  and  peculiarities  of 
the  hydrophone’s  temporal  response  cannot  be  ruled  out. 
The  model  could  be  refined  by  computing  the  reflection- 
transmission  factor  B  along  a  path  which  accounts  for  the 
finiteness  of  the  hydrophone  and  source  distances,  z  and  zp . 

Interpretation  of  the  experiments  was  simplified  be¬ 
cause  the  short  duration  of  the  incident  tone  burst  made 
possible  the  temporal  separation  of  echoes  with  different  pa¬ 
rameters  («,  m,  / )  in  most  cases.  In  true  steady-state  scatter¬ 
ing  the  axial  and  glory  echoes  overlap  in  time.  The  mutual 
interference  of  these  echoes  would  be  highly  dependent  on 
ka  when  ka>  1.  The  interference  of  glory  and  axial  rays  due 
to  several  different  parameters  is  the  plausible  cause  of  the  ka 
dependent  structures  which  are  present  in  Flax's  exact  com¬ 
putations  of  backscattering  from  large  aluminum  spheres.* 
The  exact  farfield  scattered  pressure  may  be  written  in  the 
following  form1-5  p{y,  r)  —  \p,  a/2r)/ where  the  total  form 
function  /’is  a  sum  of  the  partial-wave  amplitudes  /, 

/T*fl,y)=  £/,  (26) 

y-o 

fj  =  [2/ka)Uj(ka)3?  j{cos  y).  (27) 
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Here  r  and  y  are  the  distance  and  backscattering  angle  rela¬ 
tive  to  the  sphere’s  center,  ^y(cos  y)  is  the  Legendre  polyno¬ 
mial,  and  Uj  is  a  complex-valued  function  of  ka,j,  and  mate¬ 
rial  parameters.1'3 

When  ka  is  large,  the  axially  focused  backscattering 
amplitude  due  to  the  (n,m,l  )th  glory  ray  can  be  attributed  to 
only  some  of  the /.  Consider  the  asymptotic  expansion  of 
the  Legendre  polynomial37  for  the  case  of  j'y  1  and  y<  1 : 

^,(cos  y)  =  (y/ sin  y)mJ0[(j  +  \)y]  +0(j~m).  (28) 

The  dependence  on  y  will  be  similar  to  that  of  Eq.  (10)  pro¬ 
vided 

j~kb  -  f  (29) 

The  recognition  that  partial  waves  can  be  attributed  to  rays 
according  to  the  value  of  the  ray’s  impact  parameter  is 
known  as  the  localization  principle  of  optical  scattering  the¬ 
ory.1415,36  In  the  “physical-optics”  approximations  devel¬ 
oped  here  for  elastic  spheres  and  in  Ref.  1 3  for  fluid  spheres, 
the  dominant  ray-optical  contributions  are  for  those  rays 
having  impact  parameters  equal  to  the  radius  b  of  a  focal 
circle.  For  values  of  j  satisfying  (29),  the  phases  of  the  Uj  in 
(27)  should  be  stationary;  consequently,  there  will  be  a  signif¬ 
icant  contribution  to  /proportional  to  J0{kby).  For  values  of 
j  not  lying  close  to  kb  —  j  (for  the  b  of  any  strong  glory  ray) 
the  phases  of  the  Uj  should  be  such  that  neighboring  /  tend 
to  cancel  each  other.  A  complete  derivation  of  Eqs.  (9)  and 
( 10)  directly  from  the  partial-wave  series  (26)  would  be  com¬ 
plicated  by  the  following  caveats,  (i)  The /in  (27)  are  phase 
referenced  to  the  center  of  the  sphere  while  Eqs.  (9)  and  (10) 
are  referenced  to  C'  in  Fig.  1.  The  two  are  related  by  the 
unimodular  /-dependent  phase  factor  exp(/£c)  where  fc  is 
given  by  Fq.  (F37). 

(ii)  Equation  (9)  includes  an  additional  /-dependent 
phase  factor  exp(/<p ).  The  /  dependence  of  the  net  phase  shift, 
£c+<P—  —  ^  [2a  +  (a  —  a)(l  —  cos /)),  must  be  pre¬ 
served  in  any  asymptotic  evaluation  of  the  partial-wave  se¬ 
ries. 

Asymptotic  expressions  for  scattering  due  to  rays  trans¬ 
mitted  thru  solid  elastic  cylinders  were  obtained  by  Brill  and 
Uberall.11  Their  method  used  saddle-point  integration  to 
evaluate  the  Watson  transformation  of  the  partial  wave  se¬ 
ries  for  cylinders.  Unlike  our  results  for  spheres,  theirs  con¬ 
tain  no  ( ku)'n  factor  since  backscattering  from  cylinders  is 
not  axially  focused. 

Our  model  predicts  that  backscattered  transient  pulses 
will  be  distorted  in  the  manner  described  by  Eqs.  (22)  and 
(B3)-(B6).  Detection  of  this  distortion  was  beyond  the  scope 
of  the  present  experiment. 
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APPENDIX  A:  EFFECTS  OF  THE  REFLECTION- 
TRANSMISSION  PHASE  SHIFTS  ON  FOCAL 
PARAMETERS 

When  a  ray’s  angle  of  incidence  in  the  water  0  exceeds 
,  the  internal  rays  are  all  shear  chords.  The  phase  changes 
due  to  transmission  or  reflection  at  each  vertex  vary  continu¬ 
ously  with  6.  Consideration  of  Fig.  1  for  a  case  in  which 
0  >  0cL  will  convince  the  reader  that  outgoing  wavelet  d  'e' 
will  be  skewed  so  that  it  is  no  longer  described  by  the  geomet¬ 
ric  focal  parameters  b  and  a.  Instead,  it  is  described  by  effec¬ 
tive  parameter  be  and  a , .  The  purpose  of  this  Appendix  is  to 
illustrate  how  \b  —  be  |  <a  and  \a  —  at  |  <a  when  ka  is  large. 
The  skewing  is  only  present  for  glory  rays  having  n  =  m.  We 
simplify  the  notation  by  letting  a  subscript  n  denote  the  eva¬ 
luation  of  angles  or  derivatives  at  the  geometric  glory  condi¬ 
tion  0  =  0„my,  henceforth,  we  let  0„  denote  0 

Recall  that  s  denotes  the  distance  from  C '  at  which  a  ray 
crosses  the  exit  plane.  Figure  2  of  Ref.  13  gives  an  expanded 
view  of  the  crossing.  We  are  interested  in  rays  having  O~0„ 
so  it  is  convenient  to  use  the  expansion 

(aii 

and  to  define  the  following  dimensionless  quantities 
O'  =  a(d0/ds)„  and  0"  =  a2(d 20/ds2)„.  Evaluation  of  Eq. 
(A5)  of  Ref.  13  at  the  glory  condition  0  =  0  =  On,  gives 
0 '  =  [q  cos  0„ ) " '  where  q  =  a /{a  —  a)  is  the  spreading  fac¬ 
tor.  Inspection  of  Table  I  shows  that  0'<  1  rad.  Numerical 
evaluation  of  0 "  shows  that  it  is  also  small;  e.g.,  the  (4,  4,  1) 
ray  has  0  '~0.035  rad  and  0 "  =  —  0.007  rad. 

Recall  that  ‘P  denotes  the  total  phase  shift  due  to  each 
vertex.  It  is  convenient  to  use  the  expansion 

<P-<P„~(0-On)>P'  +  {(0-0„)2P\  (A2) 

where  —  V (s  =  b )  in  the  notation  of  Sec.  Ill, 
IP'  =  [d*P /d0)„  and  IP"  =  (d2,P/d92)„.  The  constituent 
derivatives  may  be  evaluated  at  each  vertex:  for  example, 


\  do  )„  \  d0  Jn 


(A3) 

where  IP^  is  the  shift  due  to  reflection  of  a  shear  wave  and 
IPWS  and  IPSW  are  transmission  (mode-conversion)  phase 
shifts.  These  shifts  were  computed  from  the  complex  reflec¬ 
tion  and  transmission  coefficients  of  a  plane  surface39  via  the 
algorithm  listed  in  Ref.  19.  The  IP’  and  IP'  were  obtained 
numerically. 

Let  q  denote  the  propagation  phase  delay  between  the 
entrance  and  exit  planes,  see  Eq.  (F7).  Inspection  of  (8)  gives 
V  —  c (s  —  b  )2/2a.  The  following  s  dependence  for  the 
total  phase  at  the  exit  plane  is  found  by  combining  this  result 
with  (Al)and  (A2). 


\V  -  * )  -  to,  -  IP*  ( A4) 


where  j  =  (s  —  b  ]/a,  *•'=—! P'0',  x"  —  kaa~'  +  A, 
a  =  a/a,  and  A  —  —  0'<P'  —  O'2 Ip'.  The  effective  focal- 
circle  radius  br  for  the  skewed  wave  is  the  value  of  s  for  which 
d(if  —  !P)/<k  =  0;  this  procedure  gives 


b,  -b  -X'  _  gV'e' 

a  x"  aA+ka 


(A5) 


The  aforementioned  procedure  gives  the  following  values  for 
the  (4,  4,  1)  ray:  !P'£'~  —0.85  and  A~0.il  rad.  Conse¬ 
quently,  the  shift  of  the  radius  is  negligible  when  ka  is  large, 
e.g.,  (be  —  b  )/a~  —  0.0019  for  ka  =  457.  These  results  are 
also  representative  of  those  for  rays  having  0„  closer  to  0a 
since  the  0  ‘  are  smaller  but  the  |  IP '  |  are  larger. 

The  skewed  wave  appears  to  originate  from  a  focal  cir¬ 
cle  which  is  a  distance  a,  =  k  [d  2{t}  —  IP  yds2]  ~ 1  behind  the 
exit  plane.  Application  of  (A4)  gives 

(a,,  —  aj/a  =  —aA  (aA  +  ka]~ '.  (A6) 

For  the  (4,  4,  1 )  ray,  the  shift  is  negligible  when  ka  =  457 
since  (ae  —  a)/a~  —  0.0020.  The  effective  parameters  in¬ 
fluence  the  scattering  by  replacing  a  and  b  in  Eqs.  (10)— <13), 
(17),  (20),  (24),  and  in  z„  =  z  +  a.  For  fused-silica  spheres, 
the  changes  are  negligible  unless  ka  S  100. 


APPENDIX  B:  FILTER  MODEL  OF  PULSE  DISTORTION 


Inspection  of  the  steady-state  scattering  amplitude, 
Eqs.  (10)  and  (17)  reduces  the  problem  of  the  on-axis  pulse 
distortion  to  answering  the  following  question:  What  is  the 
general  temporal  effect  of  a  filter  whose  response  to  the  func¬ 
tion  exp(  —  icoi )  is  SaSt  exp(  —  uot )?  HereS„  and  Sr  are  the 
following  frequency  responses: 

Se  (a)  =  exp  [  ie  sgn(w)  ] ,  (B 1 ) 

Sw  (a))  =  \a>  | 1  /2exp  [  -  i\ir  sgn(w)  ] ,  (B2) 

where  the  function  sgn(<w)  =  \oi\/w  and  e  is  given  by  (21).  In 
(B 1 )  and  (B2),  the  domain  has  been  extended  to  include  nega¬ 
tive  values  of  eo  by  requiring  that  phase  advancements  for 
positive  values  of  to  are  also  advancements  when  ta  becomes 
negative.  The  scattering  process  is  modeled  as  a  filter  with  a 
response  Sf  followed  by  one  with  a  response  S„.  Let  p,(t ) 
denote  the  input  to  the  first  filter;  the  output  pt  (f )  is24  23-40 


Pe(* )  =P/U  )cos  f  +pH{t  )sin  e, 
where  pH  is  the  Hilbert  transform  of  p, 


Pn(t ) 


=  1  Pf 

7T  J- 


(t-f) 


■dt' 


(B3) 

(B4) 


and  P  denotes  the  Cauchy  principal-value  integral. 

It  is  convenient  to  rewrite  the  response  of  the  second 
filter  as  Sw  —  (  —  ia>)U2  where  the  negative  real  a>  axis  is  cho¬ 
sen  to  be  the  branch  cut  in  the  evaluation  of  the  radical.  To 
obtain  the  temporal  behavior  of  this  filter  we  define  the  fol¬ 
lowing  half-order  derivative  operator2® 


(*r*)-f  pif)lirit-t')]-U2dt',  (B5) 

whereof  ’)  =  dp{t  ')/dt By  direct  evaluation  of(B5)it  may  be 
demonstrated  that 


(^)1/Je-fa"  =  Swe-‘“’.  (B6) 

Application  of  (B6)  to  each  of  the  spectral  components  of  the 
input  p(t )  to  the  second  filter,  gives  the  following  output: 

=  (B7) 

Taking/)  to  hep,  lets/>„  manifest  the  combined  effect  of  both 
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processes. 

Computations  of  the  half-order  derivative  of  a  tone 
burst,  (d  fdt  )1/2p,  with  thep,  from  (23),  show  that  the  princi¬ 
pal  distortion  is  to  the  leading  and  trailing  half-cycles. 
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I.  Introduction 


This  paper  gives  a  calculation  of  the  temporal  response  of  certain 
acoustical  systems  to  inputs  consisting  of  bursts  of  sine  waves.  The 
approximate  frequency  responses  of  the  systems  considered,  Eq.  (3)  below,  are 
characterized  by  a  45  deg.  phase  advance  together  with  a  modulus  which 
increases  with  the  angular  frequency  to  in  proportion  to  u£.  To  describe 
the  temporal  behavior  of  such  systems,  Lighthill1  used  the  following  half- 
order  derivative  operator  on  the  time  domain  function  s(t) 


r  t 


Cd/dtr  s( t) 


s(T)[(t-T)ir]“*  dT 


Cl) 


where  s(t)  =  ds/dt  (evaluated  at  t  a  t).  By  direct  evaluation  of  the 

integral,  the  steadi-state  frequency  response  h^  associated  with  this 
1  2 

operator  is  * 


(d/dt)*  e"iwt  a  h^  e’lwt, 

h„  =  <-i«0*  a  M*  e-ii^gnCco) 


(2) 

(3) 


where  the  slgnum  function  is  sgnCu  >0)  a  1,  sgnCco  <  0)  a  -1,  and  sgn(0)  a  0. 
The  sine-wave  bursts  to  be  considered  are  N  cycles  in  length  where  N  is  an 
Integer  or  half  odd-integer  so  the  s(t)  are  of  the  form 


*N(t)  a  H(uit)  H(2irN  -  cut)  sinut  (4) 

where  the  Heaviside  unit-step  funotion  H(y  <  0)  a  0,  H(y  >0)  a  1,  and 
HCy  a  0)  a  i.  Throughout  this  paper  we  take  ui  >  0. 

Section  II  reviews  several  acoustical  systems  for  which  h^  is  an 
approximate  transfer  function  for  a  range  of  u>.  This  response  is  associated 
with  the  propagation  from  two-dimensional  sources1  and  with  the  effects  of 


3 


2 

diffraction  in  axially- focused  scattering.  All  systems  considered  are  taken 
to  be  linear  so  that  signal  distortions  due  to  finite  amplitude  effects  are 
omitted.  The  effects  on  signal  shape  of  the  absorption  of  propagating  sound 
are  also  omitted.  The  examples  considered  here  constitute  one  class  of 
intrinsic  distortion.  Other  examples  of  linear  problems  for  which  distortions 
have  been  analyzed  include  those  due  to  frequency  Independent  phase  shifts^ 

4 

(caused  e.g.  by  simple  caustics  or  bottom  reflections  at  grazing  angles  of 
incidence^)  and  those  due  to  the  focusing  in  stratified  media. ® 

For  the  transients  given  by  (4).  the  integral  in  (1)  is  evaluated  in 
Sec.  III.  The  integral  also  converges  for  other  differentiable  s  for  which 
(-t)’*s(t)  ■*  0  as  t  -*■  however,  the  general  distorting  effects  of  this 
operator  will  not  be  considered  here.  Though  inputs  proportional  to  (4)  are 
not  exactly  realizable  in  physical  systems,  the  distortions  calculated  should 
be  useful  in  understanding  properties  of  real  systems.  It  is  found  here  that 
the  distortion  is  most  significant  near  the  leading  and  trailing  edges  of 
(d/dt)^s^.  The  operator  defined  in  (1)  is  a  special  case  of  the 
"semiderivative  operator"  which  has  been  used  for  solving  the  diffusion 

7  a  A 

equation  and  modeling  electrochemical  experiments.  Values  of  (d/dt)  s  for 
some  other  s(t)  can  be  inferred  from  tables  in  Sec.  7.6  of  Ref.  7  and  in 
Ref.  8.  Notice  that  when  s  =  exp(-iwt),  (d/dt)^  (d/dt)^s  *  ds/dt  so  the 
terminology  half-order  derivative  is  Justified. 

II.  Some  Acoustical  Applications  of  (d/dt)^  s(t) 

A.  Pressure  radiated  from  two-dimensional  sources 

Lighthill1  considered  the  aooustio  pressure  p(r,t)  generated  by  an 
infinite-line  thin  uniform  source  at  a  perpendicular  distance  r  from  the 
line.  The  mass  outflow  per  length  of  souroe  is  s(t),  which  has  units 
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g  cm**1  sec**1.  Integrating  the  radiation  from  each  line  element  gives  the  well 
known  result 

p(r,t)  a  (2ir)-1  j  (r^  -  r^)“^  s(t  -  r/c)dr  (5) 

■*r 


where  c  is  the  speed  of  sound  in  the  medium.  This  may  be  simplified  when  s 

-2  2 

is  a  transient  of  duration  At.  In  the  far-zone  where  r  >>  cAt,  (r  -  r  ) 
may  be  replaced  by  2r(r-  r)  so  that  (5)  becomes1 

p(r,t)  a  (c/8irr)^  (d/dt)^  s(t  -  r/c)  .  (6) 

The  condition  on  r  may  be  written  as  r  >>  X  where  X  is  the  wavelength  of 
the  principal  spectral  component  of  s(t). 

Q 

Landau  and  Lifshitz  calculate  the  radiation  from  a  radially  vibrating 
infinite  cylinder  which  has  a  cross-sectional  area  W(t).  The  unperturbed 
density  of  the  surrounding  fluid  is  PQ.  Their  result,  Eq.  (73.15),  is 

equivalent  to  (5)  with  s(t)  *  PQ  W  and  their  far-zone  result,  Eq.  (73.17), 


reduces  to  (6).  For  (6)  to  be  applicable  it  is  necessary  that 


r  »  \  »  W*. 


The  continuity  equation  for  density  may  be  used  to  show  that  the  effective 

e 

mass  outflow  per  length  is  Indeed  pQ  W. 

B.  Photoacoustic  generation  of  pulsed  sound  in  a  weak  optical  absorber 


In  a  review  of  pulsed  optoacoustic  spectroscopy,  Patel  and  Tam  apply 
the  Landau  and  Lifshitz  result  equivalent  to  (5),  to  the  calculation  of 
radiated  pressure.  The  mechanism  by  whloh  the  pressure  is  generated  is  the 
heating  of  fluid  due  to  the  absorption  of  a  pulsed  thin  light  beam  of  radius 
R.  The  electrostrictive  mechanism  is  neglected  in  their  (and  in  the 
following)  prediction  because  it  is  usually  small.11  Consider  cylindrical 
coordinates  for  which  the  light  beam  propagates  along  the  z  axis.  Let  the 
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optical  beam  power  at  z  be 


£(t,z) 


£(t,0)e' 


-az 


2tt 


[*I(t,z, 
J0 


r)rdr 


where  zl(t,z,r)  is  the  optical  Poynting  vector  at  a  radius  r  from  the 

•M 

center  of  the  assumed  symmetric  beam;  a  is  the  optical  absorption  coefficient 
of  the  fluid.  In  the  notation  of  Eqs.  (5)  and  (6),  the  Patel  and  Tam  source 
strength  prescription  becomes,  from  Eq.  (33)  and  (35)  of  Ref.  10, 


s(t)  =  a85(t,0)C“1  (7) 

P 

where  8  is  the  thermal  coefficient  of  volume  expansion  and  Cp  is  the 
specific  heat  at  constant  pressure.  For  this  prescription  to  be  applicable  it 
is  necessary  that  beam  radius  be  such  that 11 

Rd  =  (4DAt)*  <  R  «  eAt  •  X  (8) 

where  RQ  is  the  thermal  diffusion  length  for  the  pulse  duration  At  of  the 
beam  and  D  is  the  thermal  diffuslvity.  In  addition  to  (8),  it  is  necessary  that 

the  pressure  field  at  r  be  essentially  two-dimensional  so  that  a  can't  be  large  . 

r  «a-1  and  z  <<a“1;  furthermore.  At  is  assumed  to  be  sufficiently  long 
that  the  optical  transit  time  across  the  source  volume  is  «  At. 

From  the  considerations  giving  (6)  as  the  far-zone  approximation  to 
(5),  it  is  evident  that  p(r,t)  is  proportional  to  (d/dt)^  £(t  -  r/c,  0) 
where  the  beam  power  £  is  evaluated  at  the  shifted  time  which  is  r/c 
earlier  than  t.  Sinoe  £  is  intrinsically  positive,  the  calculation  of 
(d/dt)^  sN  is  only  applicable  to  photoacouatios  when  N  »  1/2. 

Readers  not  interested  in  focused  scattering  may  omit  Sec.  IIC  and  IIP 

without  loss  of  continuity. 
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C.  Axially-focused  backscattering  from  spheres 

2 

Asymptotic  approximations  to  waves  backscattered  from  large  elastic 
12 

and  fluid  spheres  have  recently  been  formulated.  The  sumary  given  here 

2 

reviews  the  connection  with  the  half-order  derivative  previously  noted.  This 
discussion  also  facilitates  the  analysis  of  the  backscattering  from  a  torus 
given  in  Sec.  IIP.  Amplitudes  due  to  glory  rays  are  enhanced  relative  to 
those  of  other  reflections  because  of  a  weak  axial  focusing.  Glory  rays  are 
defined  to  be  backscattered  rays  which  have  non-zero  impact  parameters;  Fig.  1 
illustrates  several  such  rays.  The  parameter  n  denotes  the  number  of 
internal  chords.  Figure  1  displays  as  dashed  lines,  rays  which  lie  on  either 
side  of  the  n  =  3  glory  ray.  The  wavefront  associated  with  this  glory  ray 
corresponds  to  the  curved  wavelet  d'e’.  This  wavefront  is  toroidal  and  it 
appears  to  emanate  from  a  focal  circle  formed  by  rotating  the  point  F^  about 
the  CC'  axis.  The  wave  is  characterized  by  the  focal  parameters  b,  the  radius 
of  the  circle,  and  a,  the  distance  of  the  circle  from  the  dashed  vertical 
plane . 

Let  Pj  exp(-iwt)  denote  the  incident  wave's  pressure  in  the  dashed 

plane  through  C'.  The  scattered  pressure  due  to  each  glory  ray  for  points  on 

the  backward  axis  (the  left  extension  of  the  CC'  axis)  is  computed  in  Ref.  2 

for  spheres  having  ka  >>  1  where  a  is  the  sphere's  radius,  k  *  2ff/X, 

X  *  2irc/uj,  and  c  is  the  sound  speed  in  the  surrounding  fluid  medium.  It  is 

necessary  that  the  distance  z  of  the  observation  point  from  C'  be  such 
3  —14 

that  »  X  b  .  The  scattered  pressure  of  the  nth  glory  wave  is 

P_  »  (p-  a/2z_)  exp[ik(z  «■  Jb2z ) ]G_  (9) 

n  x  n  n  n 

Gn  *  (ka)^  E  B  expCKn  ♦  e  -  iff)]  (10) 

E  *  2b(2ff<Vq)^  *”3/2  (u) 
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which  correspond  to  Eqs.  (17),  (10),  and  (13)  of  Ref.  2,  respectively; 

z  s  z  ♦  a,  B  is  the  product  of  plane-surface  reflection  and  transmission 
n 

factors  for  each  vertex  [Eq.  (7)  of  Ref.  2],  and  the  wavelet  spreading  factor 
q  is  the  limit  of  the  ratio  of  the  arc  lengths  of  d'e*  and  de  as  the 
latter  vanishes.  The  phase  n  Is  the  propagation  delay  along  the  glory  ray. 
The  phase  £  =  u  -  4f  is  independent  of  ka  (assuming  ka  »  1);  y  is  the 
combined  phase  of  the  vertex  reflection  an  transmission  coefficients.  (¥  is 
usually  0  or  ir,  the  exception  being  elastic-sphere  glories  for  the  case  in 
which  all  the  chords  are  shear  wave  and  the  angle  of  incidence  exceeds  the 
longitudinal  critical  angle.  )  The  phase  u  is  -ir/2  times  the  number  of 
internal  foci  crossed.  Three  foci  are  crossed  for  the  n  »  3  ray  and  these  are 
labeled  L1  and  L 2  in  Fig.  1. 

As  a  consequence  of  the  (ka)^  factor  in  (10),  and  pR  diverge  as 

ka  •*  ®.  This  divergence  is  due  to  the  axial  focusing.  The  geometrical  cause 

of  this  focusing  is  discussed  in  Ref.  12-14. 

The  phases  in  (9)  and  n  in  (10)  are  «  k;  they  give  rise  to  a  time 

delay  t  for  propagation.  Let  p^(t)  be  the  pressure  of  a  transient  wave 

incident  on  the  dashed  vertical  plane.  When  |e|  »  0  or  it,  the  frequency 

2 

response  in  (9)  and  (10)  is  like  that  of  (3)  so  the  transient  response  is 


Pn(t) 


3/2  „ 

m  a  EB 

’ITT 

n 


(— )^  p 
vdt;  P1 


(t  -  T) 


(12) 


where  ♦(-)  corresponds  to  lei  of  0(w).  For  intermediate  values 
of  e  there  is  an  additional  distortion  and  the  pressure  Includes  a  term 
given  by  the  Hilbert  transform*'^  of  this  pn(t). 

3/2 

The  salient  feature  of  axially-focused  backscattering  is  the  a 
proportionality  in  (12)  whereas  for  simple  specular  reflection  the 
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pressure  *  a.  With  appropriate  focal  parameters,  (9)  -  (12)  are  applicable  to 

liquid-filled  spherical  shells  as  well  as  elastic  and  liquid  spheres.  The 

internal  absorption  is  assumed  to  be  negligible.  The  cause  of  large  target 

strengths  observed  or  computed  for  certain  liquid-filled  shells1^’ ^  is  axial 

focusing.  This  is  evident  since  reported  temporal  signatures1^’ ^  show 

enhancements  (relative  to  the  specular  echo)  at  delays  of  glory  rays.  For 

example,  Fig.  7  and  8  of  Ref.  15  display  measured  signatures  and  paths  for  a 

liquid-filled  shell  having  an  acoustical  refractive  index  M  *  1.81.  The 

dominant  echo  is  from  an  n  *  2  glory  ray  which  has  2  internal  foci.  (There  is 

one  of  type  and  one  of  type  see  e.g.  Fig.  1  of  Ref.  13.) 

Consequently  u  *  -ir  and  (from  the  impedances)  !>  0  so  that  e  *  -ir  .  The 

(d/dt)^  sN  given  in  Sec.  Ill  should  describe  the  backscattering  of  sine-wave 

bursts,  however  the  data  in  Ref.  15  do  not  faciliate  comparison. 

Liquid-filled  shells  have  potential  uses  as  sonar  calibration  targets 

and,  due  to  their  large  target  strengths,  as  passive  navigational  beacons. 

For  certain  N  values,  G  may  Increase  even  more  rapidly  than  (ka£  so 

n 

(12)  is  not  applicable.  One  such  case  is  the  rainbow-enhanced  glory 

12  2/3 

predicted  to  occur  when  M  *1.18;  in  that  case  G^  *  (lea)  .  Another 

example  is  the  backward  focusing  of  the  rear  axial  reflection  which  should 

occur  when  M  *  2.  In  that  case,  elementary  diffraction  and  lens  theory 

predicts  G2  «  ka  when  ka  >>  1. 

D.  Axially-focused  backscattering  from  a  rigid  torus 

Baekscattered  pressures  proportional  to  (d/dt)^  p^  are  not  unique  to 
spheres.  They  should  also  exist  for  other  shapes.  In  this  section  I  extend 
the  treatment  in  Sec.  IIC  to  a  new  case:  backscattering  from  a  large  rigid 
immovable  torus. 


31 


9 


Figure  2  shows  the  problem  under  consideration.  The  surface  of  the 
torus  is  generated  by  rotating  the  center  of  a  circle  Cj.  about  the  point  C. 
The  radius  a  of  the  circle  is  assumed  to  be  less  than  the  distance  b  which 
separates  from  C.  In  the  analysis  which  follows  it  is  assumed  that 
lea  »  1  and  that  kb  »  1.  It  is  assumed  that  the  direction  of  propagation 
of  the  incident  plane  wave  is  -2. 

The  reflection  of  the  Incident  wave  from  the  rigid  surface  gives  rise 

to  a  virtual  backward-facing  toroidal  wavefront  shown  as  a  dashed  curve  in 

Fig.  2.  The  focal  circle  for  this  wave  (which  is  at  F^>  lies  a  distance 

a  s  a/2  behind  the  dashed  vertical  plane.  This  follows  from  geometrical 

2  12 

considerations  similar  to  those  *  which  locate  the  virtual  point-like  source 

of  the  first  reflection  (from  a  sphere)  at  AQ  in  Fig.  1.  The  radius  of  the 

focal  circle  is  b  and  the  spreading  factor  q  is  unity.  The  reflection- 

transmission  factor  B  is  unity  for  the  rigid  torus;  for  real  materials  it  is 

the  plane-surface  reflection  coefficient  for  normal  incidence. 

In  this  analysis,  multiply-reflected  glory  rays,  such  as  the  dotted 

ray  in  Fig.  2,  will  be  neglected.  Their  contribution  to  the  total 

backscatterlng  has  been  estimated,  Eq.  (4.7)  of  Ref.  17,  and  it  is  small  for 

the  large  torus  considered  here  and  especially  for  a  <<  b. 

Let  the  polar  coordinates  (relative  to  C')  of  the  observation  point 

be  r  and  y  .  A  polar  angle  of  Y  *  0°  corresponds  to  backscatterlng. 

2 

Consider  the  far  zone  in  which  r  »  kb  .  The  diffraction  integral  for  this 
single  bounce  glory  wave  is  identical  to  one  evaluated  in  Sec.  Ill  of  Ref.  12. 
When  the  pressure  incident  on  the  dashed  vertical  plane  is  Pj  exp(-lut),  the 
scattered  pressure  is 

pQ  «  (pj  b/r)frka)*  J0(u0)ei(kr  ***  "  iir)  (13) 
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where  the  argument  df  the  Bessel  function  J,,  is  :  kb  sin  y  and 

<Fa  -  ka(1  -  cosy)/2.  The  backscattering  cross  section  resulting  from  (13)  is 

equivalent  to  that  for  a  "thick-wire  loop"  in  the  corresponding 

17  18 

electromagnetic  problem.  '  The  frequency  response  of  (13)  with  Y  a  0  is 
like  that  of  (3). 

A  more  general  result  which  Includes  the  near  zone  (zJ  »  ^  b  )  is 
similar  in  form  to  Eqs.  (9)  -  (12).  The  transient  response  is  given  by  (12) 
with  the  upper  sign  and  with  the  aforementioned  focal  parameters.  These 

r  2-1 

give:  B  s  1,  E  s  2bir*/a,  zr  a  z  ♦  i  a,  and  ct  ■  a  +  i  b  zq  .  Here  z  ia 

che  distance  of  che  on-axis  observation  point  from  C' . 

III.  Half-order  Derivative  of  a  Sine-Wave  Burst 

The  half-order  derivative  of  Eq.  (4)  may  be  obtained  from 
I  (t)  =  (d/dt£  s  (t)  where  s  (t)  a  Hfot)  slnut  is  a  burst  of  infinite 

»  09  00 

duration  which  commences  at  t  a  0.  To  evaluate  ( 1 ) ,  notice  that  the  second 

term  of1^  s  a  H(ut)  u  cosut  ♦  sinut(d/dt)  H(ut)  vanishes  since 
00 

(d/dt)H(ut)  a  ufi(u)t)  where  <5  (u t)  is  a  delta  function;  consequently 

I^t)  a  J^(t-x)”^  cosux  dx  (14) 

for  t  >  0  and  I  (t<0)  a  0.  Changing  the  integration  variable  to 
y  a  (2u/ir)^  (t-x)^  gives 

I  (t)  a  (2u)^  H(ut)tC(u)  cosut  ♦  S(u)  slnut]  (15) 

OD 

C(u)  a  |U  cosC^y2)^,  S(u)  «  j"  sin(|y2)<fr,  06) 

where  u  a  (2ut/ir)^  and  C  and  S  are  Fresnel  integrals.  Define  the 

20 

following  auxiliary  Fresnel  integral  function 
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g(u)  *  Ci  -  C(u)]cos(y  u2)  ♦  Ci  -  S(u)]sin(|  a2)  (17) 

and  Eq.  (15)  becomes 

I^t)  *  H(ut)ca^[sin(wt  in)  -  2^  g(u)]  .  (18) 

It  is  well  known  that  as  u  -*■  °°,  C(u)  J,  S(u)  +  J  so  that  g(u)  -*•  0  .  It  is 
evident  from  (18)  that  for  large  cat,  I  approaches  a  sine  wave  which  is 
phase  advanced  by  45  deg  and  altered  in  magnitude  by  cj^.  These  steady-state 
properties  are  evident  directly  from  (3).  The  expressions  in  (18)  and  (15) 
are  in  agreement  with  tabulated  semiderivatives  listed  in  Ref.  7  and  8 
respectively. 

Equation  (4)  may  be  written  as  a^Ct)  *  sjt)  +  sjtt  -  <a**12TrM) .  Here 
and  below,  for  the  upper  (lower)  sign,  N  is  Integer  (half  odd-integer).  From 
the  linearity  of  the  half-order  derivative  operator,  it  follows  that 

(d/dt)*sM(t)  =  I«(t)  +  I^t  -  uf^M).  (19) 

IV.  Discussion 

For  discussing  the  Implications  of  (19),  it  is  convenient  to  introduce 
the  notations  1^ Ct )  *  (d/dt)^»N(t)  and  JN(t)  *w"^Ijj(t).  Figures  3-5 
display  calculated  JN  for  N  *  4,1,  and  i,  respectively.  The  calculations 
used  the  fast  rational  approximation  to  g(u)  given  in  Ref.  20  as  Eq. 

(7.3.33) .  The  errors  due  to  this  approximation  are  estimated  to  be  negligible 
on  the  soale  of  these  plots.  When  N  is  large,  JN  closely  resembles 
sln(ta>t  ♦  in)  except  near  its  leading  and  trailing  edges;  this  is  espeoially 
evident  on  plots  (not  shown  here)  comparing  jN  with  sin  (cat  ♦in).  The 
leading  maxima  and  minima  of  are  *  0.852,  -1.032,  0.984,  and  -1.010. 

From  the  fora  of  (19),  these  apply  to  other  for  normalized  times 
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20 

T  3  <jj t/2ir  <  N.  From  aeries  expansions  of  the  Fresnel  integrals  in  (15),  it 

can  be  shown  that  as  T  0,  so  that  Jy(t)  vanishes  with  an 

infinite  slope  at  t  *  0.  For  a  viscous  fluid,  the  leading  edge  of  the 

pressure  must  be  rounded.  For  one-dimensional  radiation  into  a  viscous  fluid, 

21 

rounding  of  the  leading  edge  has  been  modeled  for  the  case  of  a  mass 
outflow  «  s  (t). 

We  now  consider  some  implications  for  the  acoustical  systems  discussed 

in  Sec.  II.  For  scattering  responses  describable  by  (12),  the  peak-to-peak 

values  of  the  response  will  closely  approximate  the  steady-state  value  after 

the  first  cycle.  For  example,  the  peak-to-peak  values  of  j^  for  the  four 

complete  cycles  are  *  1.885,  1.994,  1.998,  and  1.999  which  are  to  be  compared 

to  the  steady-state  value  of  2.  Figures  in  Ref.  3  show  that  the  central  peak- 

to-peak  values  of  Hilbert  transformed  multicycle  sine-wave  bursts  lie  close  to 

2.  These  combined  observations  justify  the  use  of  measured  peak-to-peak 

pressures  to  confirm  quas 1-steady-state  amplitude  predictions  for  the  acoustic 

2 

glory  of  elastic  spheres. 

Figure  5  applies  to  the  photoacoustic  pressure  radiated  from  a  thin 
light  beam  for  which  the  power  $(t,0)  «  s^(t)  and  the  other  conditions  noted 
in  Sec.  I IB  hold.  The  maximum  and  minimum  of  J^(t)  we  »  0.852  and  -0.75, 
respectively. 

The  frequency  response  (3)  vanishes  at  .zero  frequency  so  the  following 

Integral 


J„(t)dt 


(20) 


must  yield  «  0  for  all  u  in  (4).  Numerical  integration  of  the 
computations  for  Fig.  3*5  gave  only  an  approximate  verification  of  this 
due  to  the  nature  of  the  algorithm  used  for  g.  It  was  desirable  to  check 
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Eq.  (19)  by  analytic  evaluation  of  A^.  For  integer  N,  the  proof  that  A^  s 
0  is  trivial  because  the  upper  sign  is  used  in  (19).  For  half-integer  values 
of  N  the  proof  requires  an  integration  of  g.  For  brevity,  we  consider  only 
the  case  N  =  i  because  the  proof  here  can  be  easily  generalized  to  R  a  3/2, 
5/2  .  .  . ;  integration  of  the  sine  terms  in  (18)  and  (19)  give 

1  OO 

u)2"*  A,  =  1-2(u  g(u)dt  (21) 

where  u  =  (2wt/ir)^.  The  integral  in  (21)  may  be  evaluated  by  first  using 
properties20  of  the  Fresnel  integrals  to  put  (17)  in  the  following  form 

g(u)  a  1“  e"iruy  sin(|  y2)dy  ■.  (22) 

Inserting  (22)  into  (21)  and  changing  the  order  of  the  integrations  yields 

22 

Aj  s  0  where  use  Is  made  of  tabulated  Integrals.  Equation  (22)  and  the 
requirement  that  the  integral  in  (21 )  yield  Aj  *  o  are  properties  of  the 
auxiliary  Fresnel-integral  function  g  which  appear  to  have  not  previously 
been  noted. 

The  part  of  for  T  >  N  may  be  referred  to  as  the  "tail"  or 

"wake"  of  the  response.  The  result  »  0  requires  the  tail  to  be  more 

9 

prominent  for  N  =  1/2,  3/2,  5/2  .  .  .  than  for  integer  N.  It  is  well  known 

that  pressure  transients  from  two-dimensional  sources  vanish  slowly  as  t  ■* 

Consider  now  the  trailing  edge  for  some  related  radiation  problems.  For 

finite- length  cylindrical  sources  with  mass  outflows  of  finite  duration,  the 

radiated  pressure  will  be  of  finite  duration  for  a  dissipationless  fluid.  J 

Radiation  from  a  piston  driven  with  a  velocity  «  s^(t)  has  been  modeled  by 
24 

Beaver  with  the  result  that  the  pressure  is  of  finite  duration.  The  long 
tall  on  should  be  suppressed  in  real  axlally-fooused  scattering  because 
(12)  is  only  applicable  for  speotral  components  having  ka  »  1. 
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Figure  Captions 

Fig.  1.  Backscattered  rays  from  a  sphere  for  which  the  acoustical  refractive 
index  M  =  0.6.  For  the  3  chord  glory  ray,  the  associated  back-facing 
wavefront  d'e'  is  shown  together  with  the  locations  (L^  and  Lg)  of  internal 
foci.  This  wave  is  axially  focused  along  the  left  extension  of  the  CC'  axis. 

Fig.  2.  Cross-sectional  view  of  the  torus  considered  in  Sec.  IID.  The  figure 
may  be  rotated  about  the  CC'  axis.  The  reflected  wave  is  axially  focused 

along  the  left  extension  of  the  CC'  axis  when  the  incident  wave  propagates  in 

the  -z  direction.  Scattering  due  to  the  dotted  ray  is  omitted. 

Fig.  3.  The  solid  curve  is  the  calculated  JN  3  u  (d/dt)7sN  for  the  4  cycle 
burst  given  by  the  dashed  curve.  The  normalized  time  is  T  *  a>t/2ir. 

Fig.  4.  Like  Fig.  3  but  for  N  =  1. 

Fig.  5.  Like  Fig.  3  but  for  N  »  1/2.  For  real  photoacoustic  sources,  the 

pressure  transitions  at  T  •  0  and  0.5  will  be  smoother  than  due  to  the 
convolution  over  the  finite  optical  beam  area  discussed  in  Ref.  11. 
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Abstract 

Expressions  are  derived  from  Hie  theory  for  the  cross-polarized  and  non-cross- 
polarized  measurement  configurations  useful  for  the  study  of  the  optical 
glories  of  transparent  spheres. 
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1  2  3—5 

Optical  backscatterlng  from  cloud  droplets,  '  bubbles  in  liquids, 

6 

and  other  transparent  spheres  is  enhanced  due  to  a  weak  focusing  along  the 
backward  axis.  This  enhanced  scattering  is  referred  to  as  the  glory.  Axial 

7 

focusing  should  also  enhance  the  forward  scattering  from  drops  and  from 

5 

bubbles.  The  purpose  of  this  letter  is  to  relate  certain  observable 

1  8 

intensities  with  Mie  theory  *  for  the  case  of  plane-polarized  incident  light. 

A  diagram  of  the  type  of  experiment  under  consideration  is  shown  in  Fig.  1.  A 

sphere  is  illuminated  with  a  beam  propagating  in  the  z  direction  having  the 

incident  E  field  along  the  x  axis.  The  near  forward  or  backward  scattering  is 

viewed  via  ideal  polarizers  which  may  be  oriented  to  either  block  or  pass  x- 

polarlzed  fields.  The  polarizers  lie  in  planes  which  are  normal  to  the  z  axis 

in  the  far-scattering  zone.  Previous  descriptions  of  the  transmitted 

intensities  based  on  Mie  theory  (e.g.  Eq.  6  of  Ref-  4)  were  nonuniform  in  that 

they  were  applicable  only  near  the  forward  or  backward  directions.  The 

description  presented  here  is  not  restricted  by  choice  of  angle. 

It  is  necessary  to  review  certain  results  of  Mie  theory.  The  coordinate 

system  shown  in  Fig.  2(a)  is  equivalent  to  the  one  used  in  Ref.  1.  At  an 

observation  point  Q  the  orthogonal  unit  vectors  »p, jig,  and  are  directed 

according  to  the  displacement  of  Q  with  infinitesimal  increases  in  r,  9  , 

and  Cf ,  respectively.  Let  a  denote  the  radius  of  the  sphere  and 

k  *  2it/Xo,  where  XQ  denotes  the  wavelength  in  the  outer  media  (which  may  be  a 

2 

liquid  in  the  case  of  scattering  from  bubbles).  In  the  far  zone  (r  »  ka  ), 
without  a  polarizer  the  scattered  field  is  <Efl  *  EqJq  ♦  E^gy  where 

Eg  *  f  S2  cos?,  *  -  f  S,  sin(P,  (1) 

f  a  (-lE^/kr)  exp(lut  -  ikr),  and  and  S2  are  the  standard  Mie  amplitude 
1  9 

functions  *  which  are  to  be  computed  for  the  appropriate  scattering  angle  0, 
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size  parameter  ka,  and  relative  refractive  index  m  of  the  sphere ,  and  E^  is 
the  incident  wave's  amplitude. 

The  polarizers  (e.g.  sheets  of  idealized  Polaroid)  do  not  redirect  the 
scattered  waves;  consequently  the  local  transmitted  field  Ej  e^  is  orthogonal 
to  £r.  The  values  for  the  parameter  J  are  defined  as  follows.  Consider  the 
effect  of  polarizer  1  on  a  z-dlrected  wave.  If  it  is  oriented  to  perfectly 
transmit  E  fields  along  e  while  perfectly  absorbing  those  along  e  .  then 
J  *  5.  If  it  transmits  E  fields  along  while  absorbing  those  along  £x, 
then  J  a  6.  (This  convention  for  j  is  Justified  below.)  These  definitions 
apply  also  to  the  effects  of  polarizer  2  on  a  wave  directed  along  the  -z  axis. 
The  unit  vector  e_  may  be  constructed  by  removing  from  e„  that  part  of  e_ 
which  lies  along  a  .  It  is 


*5 


*  “sir*  *  VC,  * 

’  »5  -  sin0  cos^f^),  (2) 


where  h_  *  |e  X*!'1  a  [1  -  (sin0  cos^)2]”1/2  and  the  right  side  of  (2) 
follow  from  expressions10  for  J»x  in  terms  of  the  «rt  Jfe  • Likewise,^  may 
be  constructed  by  removing  from  e„  that  part  of  e  lying  along  e  .  It  is 

^  y  fmr 


it  •  "6  St  *  <Ar  Xir>  *  V*  -  (3) 

_i  2  —1/2 

where  h,  »  |e  x  el  *  Cl  -  (sln9sinf)  ]  .  In  Fig.  (2a),  e*  is  shots), 

o  n»y  a*r  *•*> 

Note  that  for  0  a  0  or  9  a  w  that  e_  a  e  and  e.  a  e. 

mZ  *0  my 

For  eaoh  case  the  propagating  transmitted  field  becomes 

E.  a  e>«  E  .  It  is  convenient  to  define  amplitude  functions  S.  such  that 
J  MtJ  J 

ej  ■  f  sj  •  J  *  <*> 

where  the  ohoioe  of  values  for  J  were  seleoted  to  avoid  confusion  with  the 
unrelated  amplitude  functions1  s3  and  3,  (which  vanish  for  spheres).  It  is 
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convenient  to  express  the  in  terms  of  the  standard  amplitude  functions 
and  Sg  via  Eqs.  (1)  -  (3).  This  procedure  gives 

Sg  s  h^ (3-2030  cos^  ♦  S1  sinfy  (5) 

Sg  s  hg  sin^P  cos^CS 5  cos0  -  S^).  (6) 

In  agreement  with  elementary  considerations,  =  ±ir/2  gives  =  S^; 

also  ^  =  0  or  it  give  s  S2  in  the  forward  hemisphere  and  Sg  s  -  S2  in 

the  backward  hemisphere. 

The  amplitudes  Sg  and  are  associated  with  "cross-polarized**  and 

"non-cross-polarized"  scattering  measurements ,  respectively. ^  (This  definition 

of  cross-polarization  is  not  universally  held.)  When  describing  their 

associated  intensities,  it  is  convenient  to  use  a  normalization  factor  Xr  3 
2  2 

I. a  /4r  where  I.  is  the  incident  intensity  and  I  is  the  Intensity  at  a 
11  r 

radius  r  »  a  for  a  perfectly  reflecting  sphere  of  radius  a  (as  predicted  by 
naive  ray  optics).  The  Intensity  transmitted  by  a  j-orlented  polarizer  is 

li  *  lr  4»sj  (7) 

For  observations  made  in  the  backward  hemisphere  with  the  apparatus  illustrated 
in  Fig.  1,  the  (partial  transmission)  losses  due  to  the  beam  splitter  reduce  the 
fields  from  those  given  by  (1).  The  final  result  (7)  is  still  applicable  if  I 

r 

is  Interpreted  to  be  the  perfect  reflector  result  with  losses  and  the 
polarization  dependence  of  the  losses  are  negligible.  K  similar  interpretation 
facilitates  the  application  of  (7)  to  scatterers  in  liquids  viewed  via  windows 
with  detectors  (e.g.  cameras)  plaoed  in  air.  In  the  latter  case,  the  6  in  (5) 
and  (6)  becomes  the  refracted  scattering  angle  while  the  S1  and  S2  are 
evaluated  for  a  scattering  angle  within  the  liquid. 
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Equations  (5)  -  (7)  are  the  principal  results  of  this  analysis.  The 
discussions  which  follow  concern  their  application  to  backward  and  forward 
scattering  from  large  transparent  spheres. 

To  facilitate  discussion  of  the  backscattering  it  is  convenient  to  use 
the  angles  illustrated  in  Fig.  2(b);  S1  'and  S2  are  evaluated  at 
6  s  ir  -  y .  Equations  (5)  and  (6)  may  be  written  as 

S5  s  JhgCs”  -  S+  cos2<f  ♦  (1  -  cosY)2S2  cos^],  (8) 

Sg  s  -  Jhg  sin2<P  CS'*’  -  (1  -  cosy)S23,  (9) 

where  S*  =  —  S2  appear  naturally  in  the  computation  of  Nle  scattering, ^ 

5  -1/9  •  9—1/9 

hj  *  [1  -  (siny  co aftl  ,  and  hg  *  [1  -  (siny  sin<p)  ]  .  The  previous 

q 

approximations  for  polarized  and  cross-polarized  Hie  scattering  are  given  by 

replacing  the  h^  by  unity  and  omitting  the  terma  proportional  to  (1  -  cosy). 

Backward  axial  focusing  is  only  significant  for  quite  small  values  of  y.  The 

angles  of  interest  have  y  <  5*  for  which  (1  -  cosy)  <  0.004  and  the 

approximations  Introduce  negligible  errors.  Evidently  the  errors  should  be 

largest  for  the  cross-polarized  case  as  y  ♦  0,  since  for  spheres1  S*  -*•  0 

though  S2^  0.  [Representative  plots  of  Ig  for  bubbles  (in  which  S2  was 

omitted  from  (9))  are  given  in  Ref.  4  and  5;  the  first  maximum  of  |s+|was  found 

to  be  near  y  of  X  /a  rad.  which  is  «  5*  in  cases  of  strong  bsckscsttering. 

o 

In  those  papers  the  azimuthal  angle  is  denoted  by  ?.]  The  exact  and 

approximate  Ig  vanish  for  y  *  0  and  for  ^  *  0,  ±ir/2,  and  it,  these 

symantrles  are  evident  in  photographa  of  cross-polarized  glory  from  bubbles 

4  ♦ 

Ig  manifests  the  glory  beoauae  that  part  of  S  due  to  unfocused  reflections 
is  negligible  when  y  <  5*. 

For  consideration  of  forward  scattering  it  is  convenient  to  rewrite  (5) 
and  (6)  as 
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Sg  =  ihgtS*  -  S'  cos2?-  (1  -  cos0)2S2  cos^p],  (10) 

S6  3  "  ihg  sin2^  [S'  ♦  (1  -  eos0)S2],  (11) 

which  show  that  the  dependencies  on  S*  and  S~  are  reversed  froa  those  in  (8) 
and  (9).  When  0  *  0,  Ig  vanishes  since1  S1  *  S2  and  S'  »  0.  For  spheres 
large  enough  for  axial  focusing  to  be  significant  (ka  >  200),  1^  and  Ig  are 
of  similar  magnitude  in  near  backscattering  but  not  in  near  forward  scattering. 
The  diffraction  contribution  to  near  forward  scattering  is  (Ref.  1,  p.  210) 

(S,)d  *  (S2)d  »  W>d  »(ka)2a'1  J^a) 

where  m  >  ka0,  For  comparison,  the  leading  contribution  to  S~  [which  from 
(11)  is  evidently  the  intrinsic  cross-polarized  scattering]  must  Increase  less 
rapidly  with  increasing  ka;  for  example,  transmitted  toroidal  wavefronts1’^ 
contribute  to  3“  (and  to  S>)  in  proportion  to  (ka)^/2.  Due  to  the  relatively 
large  magnitude  of  (S2)d,  errors  due  to  omitting  the  term  proportional  to  S2 
should  be  more  significant  in  (11)  than  in  (9).  Forward-directed  toroidal 
wavefronts  should  manifest  themselves  in  measured  Ifi;  due  to  the  non-ideal 
aspects  of  polarizers  and  sources,  care  should  be  taken  to  remove  (via  a  spatial 
filter)  the  strong  unscattered  beam. 

This  work  was  supported  by  the  Office  of  Maval  Research  and  by  an  Alfred 
P.  Sloan  Foundation  Fellowship. 


49 


references 


6 


1.  H.  C.  vaa  de  HulsC,  Light  Scattering  by  Small  Particles  (Wiley,  New  York, 
1957). 

2.  V.  Khare  and  H.  M.  Nusaenzvelg,  "Theory  of  the  glory,"  Fhys.  Rev.  Lett. 

38,  1279-1282  (1977). 

3.  D.  S.  Langley  and  P.  L.  Mars  ton,  "Glory  in  optical  backscattering  from  air 
bubbles,"  Phys.  Rev.  Lett.  47,  913-916  (1981). 

4.  P.  L.  Marston  and  D.  S.  Langley,  "Glory  in  backscattering:  Mie  and  model 
predictions  for  bubbles  and  conditions  on  refractive  index  in  drops," 

J.  Opt.  Soc.  Am.  ]2_,  456-459  (1982);  In  the  final  section,  the  equation 
pT  ■  -1  should  be  px  ■  1. 

5.  P.  L.  Marston,  D.  S.  Langley,  and  D.  L.  Kingsbury,  "Light  scattering  by 
bubbles  in  liquids:  Mie  theory,  physical-optics  approximations,  and 
experiments,"  Appl.  Sci.  Res.  38,  373-383  (1982);  P.  L.  Marston,  "Light 
scattering  by  bubbles  in  liquids:  comments  and  application  of  results  to 
circularly  polarized  incident  light,"  Appl.  Sci.  Res.  40,  3-5  (1983). 

6.  J.  J.  Stephens,  P.  S.  Ray,  and  T.  W.  Kitterman,  "Far- field  Impulse 
response  verification  of  selected  high-frequency  optics  backscattering 
analogs,"  Appl.  Opt.  14,  2169-2176  (1975). 

7.  H.  M.  Nussenzveig  and  W.  J.  Wiscombe,  "Forward  optical  glory,"  Opt. 

Lett.  5,  1279-1282  (1980). 

8.  G.  Mie,  "Beit rage  zur  Optik  t ruber  Median,  spezlell  kolloidaler 
Metallosungen,"  Ann.  Phys.  (Leipzig)  25,  377-445  (1908).  [English 
Translation  No.  79-21946  (National  Translation  Center,  Chicago,  1979)]. 

9.  W.  J.  Wiscombe,  "Improved  Mie  scattering  algorithms,"  Appl.  Opt.  19. 
1505-1509  (1980). 

10.  G.  Arfken,  Mathematical  Methods  for  Physicists.  2nd  ed.  (Academic,  New 
York,  1970),  p.  84. 


50 


Fig.  2.  Spherical  coordinates  used  in  the  description  of  scattering  to  a 
point  Q  in  the  forward  hemisphere  (a)  and  backward  hemisphere  (b).  The 
sphere  is  centered  on  0.  Tlte  polarizers  are  in  planes  perpendicular  to  the 
t  z  axes  such  that  they  interrupt  the  lines  0Q  when  r  is  large. 
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